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ABSTRACT 

We present interferometric angular diameter measurements of 21 low-mass, K- 
and M- dwarfs made with the CHARA Array. This sample is enhanced by adding 
a collection of radii measurements published in the literature to form a total data 
set of 33 K-M dwarfs with diameters measured to better than 5%. We use these 
data in combination with the Hipparcos parallax and new measurements of the 
star's bolometric flux to compute absolute luminosities, linear radii, and effective 
temperatures for the stars. We develop empirical relations for ~K0 to M4 main- 
sequence stars that link the stellar temperature, radius, and luminosity to the 
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observed {B -V), {V - R), {V - 1), {V - J), {V - H), {V-K) broad-band color 
index and stellar metallicity [Fe/H]. These relations are valid for metallicities 
ranging from [Fe/H] = —0.5 to +0.1 dex, and are accurate to ~2%, ~5%, and 
~4% for temperature, radius, and luminosity, respectively. Our results show that 
it is necessary to use metallicity dependent transformations in order to properly 
convert colors into stellar temperatures, radii, and luminosities. Alternatively, we 
find no sensitivity to metallicity on relations we construct to the global properties 
of a star omitting color information e.g., temperature-radius and temperature- 
luminosity. Thus, we are able to empirically quantify to what order the star's 
observed color index is impacted by the stellar iron abundance. In addition to 
the empirical relations, we also provide a representative look-up table via stellar 
spectral classifications using this collection of data. Robust examinations of sin- 
gle star temperatures and radii compared to evolutionary model predictions on 
the luminosity - temperature and luminosity - radius planes reveals that models 
overestimate the temperatures of stars with surface temperatures < 5000 K by 
~ 3%, and underestimate the radii of stars with radii < 0.7 R© by ~ 5%. These 
conclusions additionally suggest that the models over account for the effects that 
the stellar metallicity may have on the astrophysical properties of an object. By 
comparing the interferometrically measured radii for the single star population 
to those of eclipsing binaries, we find that for a given mass, single and binary 
star radii are indistinguishable. However, we also find that for a given radius, 
the literature temperatures for binary stars are systematically lower compared to 
our interferometrically derived temperatures of single stars by ~ 200 to 300 K. 
The nature of this offset is dependent on the validation of binary star temper- 
atures; where bringing all measurements to a uniform and correctly calibrated 
temperature scale is needed to identify any influence stellar activity may have on 
the physical properties of a star. Lastly, we present a empirically determined HR 
di agram using fundamen tal properties presented here in combination with those 



m 



Boyajian et al.l (120121 ) for a total of 74 nearby, main- sequence. A- to M-type 



stars, and define regions of habitability for the potential existence of sub-stellar 
mass companions in each system. 



Subject headings: Stars: fundamental parameters. Stars: late-type. Stars: low- 
mass. Infrared: stars. Techniques: interfero metric. Techniques: high angular 
resolution. Stars: atmospheres. Stars: general, (Stars:) Hertzsprung- Russell and 
C-M diagrams, (Stars): planetary systems 
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Introduction 



Direct size measurements of low-mass stars represent vital tests of theoretical rn o dels o f 
stellar evolution, structure, and atmospheres. As seen in the results of iBerger et al.l ( 120061 ) . 
notable disagreements exist between interferomet rically determined r a dii an d t hose calcu- 



l ated in low-mass stellar models such as those of IChabrier &: Baraffd ( 119971 ) or ISiess et al. 



( 119971 ) in the sense that interferometrically obtained values for the stellar diameters are sys- 
tematically larger by more than 10% than those predicted from models. A similar trend 
exists for the linear radii and temperatures of low-mass stars obtained from eclipsing bi- 
naries (^EBs), wh ich are systematically larger and cooler (respectively) for given mass (see 
Irwin et al.l |20 111 , and references t herein). Additional evidence fo r this ongoing discrepanc y 



can be found in the work s of, e.g., iLopez-Morales fc RibasI (|2005[ ). Ivon Braun et al.l ( 120081 ) 



and lBoyajian et al.l (|2008[). providing motivation for adjus tments to models in order to match 
the observations (IChabrier et al.ll2007l : iDotter et al.ll2008l ) and fulfilling the obvious need for 
more directly determined diameter s of late-type dwarfs to provide further constraints to 
stellar models ( iBoyajian et al.ll2010( ). 



The emergence of the observed radius discrepancy over the years (e.g. iPopperl Il997l : 
Torres fc RibasI |2002| ) for low- mass EBs has triggered different theoretical scenarios to rec- 
oncile model predictions with observations. iRibas et al.l ( l2008l ) identify that the elevated 
activity levels in b inary stars mu s t be a ccounted for in defining the properties of low-mass 
stars. In addition, iMorales et al.l (120 lOf ) show that observational analysis, and thus the de- 
rived parameters of a system, may be biased if the activity (manifested as star spots) is not 
properly accounted for in EB light curves. Strong magnetic fields wit hin the stellar interior 
have been shown to inhibit convection in fully convective M-type stars (jMuUan fc MacDonald 
200ll ). Other contributing factors to the EB radii offsets with respect to models such as the 



orbital period an d the rotation rate have been considered and refuted (IKraus et al.l 12011 
Irwin et al.lboill ). 



While such proposed scenarios alleviate the discrepancy between models and observa- 
tions of binary stars, they are not capable of explaining the presence of the same shortcomings 
observed in predictions of sing le star properties. However, for the case of single stars, the 
interferometric observations in iBerger et al.l (|2006[ ) suggest that theoretical models for low- 
mass stars may be missing some opacity source in or der to properly repro duce observed 
radii. This result was challenged by the conclusions in iDemory et al.l ( 120091 ). which reveal 
no such correspondence. In many instances in the literature, the effective temperatures of 
stars have been shown to be over estimated by models, also alluding to a possible connection 



with stellar abundances (e.g., see lBoyajian et al.ll2008l and references therein). 



The first measurement of a stellar diameter for a star other than our Sun was that of 
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Betelgeuse in 1921 ( iMichelson &: Peaselll92l[ ). Since then, several hundred stars have been 
observed at very high spatial resolution to measure their angular sizes. The most difficult to 
target are stars on the low-mass end of the main sequence. Although plentiful and nearby, 
these stars are intrinsically small and faint, and require sensitive high angular resolution 
techniques to resolve them adequatel}0. The current number of single K- and M- dwarfs 



whose angular s izes have been ob t ained with long- 



preci s ion is 17 ([Lane et al. 



boOll: ISegransan et al. 



jaseline interferometry to b e tter than 5% 

2003 iBerger et al. 



2003 



Di Folco et al. 



2006 : iBovaiian et al. 



20091: 



200 



von Braun et al. 



Kervella et al. 



2011 



20081: 



Demory et al. 



2009 



Ivan Belle fc von Braun 



2OI2I ) and fewer than half of these are M-dwarfs. Empirical 



radius and temperature measurements toward the low-mass end of the main-sequence are 
sparsely available, and any placement of low-mass stars onto an HR diagram presently relies 
upon model atmosphere results that are in need of observational verification and constraints. 
We also identify that with the limited amount of the data available in the literature presently, 
it is difficult to discern whether or not the discrepancy between models and observations can 
be explained by a unified theory to encompass issues observed in both binary and single 
stars. 

In the context of the present paper, we stress the observational components of an ex- 
tensive interferometric survey to determine the fundamental properties of ne arby stars. The 



results of our recent paper focus on main-sequence A-, F-, and G-type stars (IBoyajian et al. 



20121 ). Here, we present high-precision interferometric observations to determine their an- 
gular diameters of 21 late-type K- and M-type dwarfs (§ In order to characterize these 
systems, we use trigonometric parallax values and literature photometry to calculate physical 
diameters and stellar bolometric fluxes. This allows us to determine their effective temper- 
atures and luminosities (§ [3]). In § HJ we take our stellar parameters and literature data to 
establish relations between stellar astrophysical parameters and observables. Throughout 
the discussion, we elaborate on the agreement of the data to predictions of stellar atmo- 
sphere and evolutionary models. We address the aforementioned radius and temperature 
discrepancy with models in § [51 We provide motivation to validate estimates of binary star 
temperatures in § [6l as they are found to be several hundred Kelvin lower than those of 
single stars, and we conclude in § [3 



^Their linear sizes amount to only several Jupiter radii, and several orders of magnitude less luminous 
than the Sun. In comparison, Betelgeuse is ^ 11, 000 times the size of Jupiter and ^ 140, 000 the luminosity 
of our Sun! 
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2. Observations with the CHARA Array: Angular Diameters 

We acquired interferometric observations at the CHARA Array with the Classic beam 



comb iner in single-baseline mode in the near-infrared K' and H bands (jten Brummelaar et al 



20051 ) for 21 K- and M-type dwarfs (Table [T]). Targets were selected to be of northern dec- 
lination ( — 10° < dec < 90°), to have no known companion with separation < 2", to be 
bright enough for the instrument (if, K' < 6.5, V < 11), and to have estimated angular sizes 
large enough to resolve the stellar disk to better than 5% precision [9 > 0.4 mas for our 
configuration). The majority of the data were collected from 2008 through 2011 with the 
longest of CHARA's baselines -Bxeiescope #i/Teiescope #2- 5si/Ei = 331 m, -Bei/wi = 313 m, 
5s2/Ei = 302 m, 5si/E2 = 279 m, 5si/wi = 279 m, 5s2/e2 = 248 m, and -Bei/w2 = 218 m. 
Additionally, we made use of CHARA archived data for GJ 687 taken in 2004, and for GJ 526 
from 2006. 

We employ the standard observing strategy for interferometric observations where each 
science target is obseryed in bracketed sequen ces with calibrator stars (for a detailed discus- 



sion see 



Boyajianll2009l : iBoyajian et al.ll2012l ). In short, we were cautious that every science 
target was observed with (1) at least two calibrators, each chosen to be unresolved point 
sources in close proximity to the object on the sky, (2) on a minimum of two nights, and 
(3) with more than one baseline configuratioii^. We collected an average of 22 bracketed 
observations (ranging from n = 2 to 51) for each star, ideal for reducing noise and increasing 
the reliability of the final diameter fits. The low number of observations for GJ 702B {n = 2) 
is a consequence of technical difficulties we had to acquire and lock tip-tilt on the secondary 
star which lies ~ 4.5 arcseconds from GJ 702A, the much brighter primary in the visual 
binary syste ng. However, the observations of GJ 702B still meet the aspects 1-3 outlined 
above. In order to estimate the calibrator star's angular diameter, 6'sed5 we fit fiux calibrated 



photo metry to a Kurucz model spectral energy distribution (for details, see iBoyajian et al. 



20121 ). As noted in item 1 above, an ideal calibrator is an unresolved point source in order to 



minimize any biases that may arise through the assumed cali brator star diameter. Simila r 



to the strategy employed in Paper I, we follow the example in Ivan Belle fc van Belld ( l2005l ). 
aiming to select calibrators with 6'sed < 0.45 mas. In this way, the errors on the object's 
calibrated visibilities are mainly a product of measurement error, and the relative errors 
introduced by calibrator angular size prediction error are minimized. All targets observed in 
this work have at least one calibrator that meets this criteria. In several cases however, we 



^The exception to item #3 is GJ 33, only observed on Sl/El 

•^The interferometric field of view is on the order of 2 arcseconds, and thus the individual component 
visibilities are not affected by incoherent light on the detector 
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use partially resolved calibrator stars with angular diameters larger than 6'sed = 0.45 mas. 
All calibrated points are cross-checked for consistency with respect to one another in order 
to identify any biases that stem from using partially resolved calibration sources, as well as 
any signs of duplicity (i.e., if a selected calibrator star is a previously unidentified binary 
star.) A log of our observations can be found in Table [1], and a list of our calibrators is 
presented in Table [2l 

To solve for the angular diar neter of our science star, we used a non-linear least-squares 
fitting routine in IDL (MPFIT, iMarkwardtl 120091 ) to fit the calibrated data points to the 
uniform disk and limb dark ened visi bility function for a sin gle star using the equations in 
Hanbury Brown et al.l (Il974l ) (see also Ivon Braun et al.ll201ll ) . To correct for limb dar kening. 



we use the linear limb- darkening coefficients from PHOENIX models presented in I Claret 



(I2OOOI ). where A represents the wavelength of observation {H or K, where here we assume 
K ^ K'). The determination of a limb-darkened angular diameter is very forgiving in the 
infrared, where the change from 6'ud to 6'ld is typically on the order of 2 — 3 %. In addi- 
tion to a wavelength dependence, the coefficients are also dependent on the star's effective 
temperature, surface gravity, and metallicity. 

When adopting limb- darkening coefficients, we find that modest errors on the assumed 
stellar parameters of temperature and gravity (ATeff = 200 K, A log (7 = 0.5 dex) contribute 
to under a tenth of a percent to the overall error budget of the limb-darkened angular diam- 
eter. The dependence on the coefficients with respect to metallicity is negligible throughout 
the range of metallicity values in our sample, so for all objects, the respective limb-darkening 
coefficients assume solar metallicity. In order to assign the appropriate limb- darkening coef- 
ficient to an object, we perform an iterati ve procedu re, beginning with an initial guess of the 
temperature based on spectral type from lCoxl (I2OOOI ). F or all stars , we adopt a surface grav- 
ity value for lower main-sequence stars to be logg = 4.5 (jCoxll2000l ). The resulting fit for the 
limb-darkened angular diameter using this coefficient is then combined with the measured 
bolometric flux (see Section 13. ip to derive a temperature. This temperature is then used to 
identify a new limb darkening c oefficie i it, if n eeded. All temperatures are simply rounded to 
the nearest temperature in the IClaretl (I2OOOI ) grid, which is in 200 K increments. 



Our results and associated uncertainties are given in Table [3|, and Figures [H [21 [31 and 
[H display the fits to the data. In summary, we measured the angular diameters of 9 K stars 
and 12 M stars with a precision of better than 2.5%. 

We intentionally chose several target stars that have previously been observed by other 
interferometers and beam combiners to check whether our measurements are consistent with 
other published values and/or to improve on previous measurements with larger uncertain- 
ties. Table [H shows this comparison for 12 of the stars presented here along with spectral 
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type, angular diameter value, and instrument used. The table further provides the relative 
error of each measurement as well as the statistical difference between our measurement 
and the previously published value. We indicate by marking a measurement with a c where 
the previously published diameter error is greater than 5%, leaving only 5 stars (7 litera- 
ture measurements) that have <5% uncertainty for comparison with our results (Figure |5]). 
Literature measurements with large errors are consistent with CHARA measurements of di- 
ameters < 1 mas. This is a direct consequence of the high angular resolution available with 
CHARA's long basehnes. 



An exception is the work in lBerger et al.l (120061 ) . where they also use the CHARA Array 
to r neasure angu l ar dia meters of stars., although their errors are still quite large. However, 
the iBerger et al.l (120061 ) results are nonetheless different from our work, possibly from their 
measurement errors being underestimated. For instance, as listed above, our observing 
strategy is to use more than one calibrator, more than one baseline, and obs erve on more 



than one night in order to minimize calibration issues that might arise. In iBerger et al. 



( I2OO6I ) this is not the case for GJ 15A, GJ 687, or GJ 880, and although GJ 526 is observed 
on more than one occasion, it lacks the additional calibrator and baselines. We also make use 
of the H-hand filter for our observations, improving the interferometers resolution by a factor 
of ~ 1.5. Upgrades to the instrument have moreover impr oved sensitivity and m inimized 
calibration effects arising in data acquisition and reduction (jSturmann et al.ll2010l ). 



ometer (PTI) (ILane et al 



Interferometer (VLTI) (jSegransan et al. 



2001 



The other measurements compared in Table HI come frora the Palomar Testbed Interfer- 

and the Very Large Telescope 



van Bel 



e fc von Braun 



200^ 



2009 



^; Demory et al. 



2OO9I ) with much shorter base- 



lines than CHARA. We see no statistical difference between the measurements at PTI, VLTI 
and our own, and because of CHARA's larger baselines, we improve the precision of these 
measurements by no less than a factor of three for all previous results, and in one case by 
up to a factor of 10. Note that iDemory et al.l (120091 ) state that beam combiners such as 
CLASSIC that do not use spatial filtering are prone to systematic calibration errors. While 
it is true that beam combiners without spatial filtering do indeed have larger errors, we 
have found no evidence either in our own measurements or in the literature that these errors 
are systematic. Furthermore, we point out that for many years no w, the CLASSIC beam 



comb iner has included spatial filtering optimized for the faint targets (jten Brummelaar et al. 
2008h . 



In conclusion, we present new interferometric diameter measurements of 21 low-mass 
stars with an average precision of 1%. Nine of these (five M-stars and four K-stars) have 
not been measured before this work. Of the 12 previously observed stars, our new data 
improve the precision for seven of these objects to better than 5%. The five remaining stars 



- 8 - 



in our sample (GJ 15A, GJ 166A, GJ 380, GJ 411, and GJ 699) are used as consistency 
checks across different interferometers and instruments. We find that the average ratio 
^LD.this work/^LD, reference = 1-008 proves that there is excellent agreement among the different 
data sets. 



3. Stellar Parameters 

In this section, we use our angular diameters in combination with literature data to 
determine stellar parameters for the 21 stars in Table [31 We use these same methods and 
procedures to determine properties of stars in the literature with interferometrically deter- 
mined radii (see Section [373]) . 



3.1. Stellar Diameters, Effective Temperatures, and Luminosities 



All of the observed stars are bright en ough to have precise trigonometric parallaxes 
measured by Hipparcos (Ivan LeeuwenI 120071 ) . These values are converted to distances d, 
which range from 1.82 to 9.75 parsecs, and are known with an average precision of 0.5%. 
The physical stellar radii R are trivially calculated via 



0^rt = 



2R 
J 



We produce stellar energy distribution (SED) fits for all targets using fiux calibrated 
photometry pubhsh e d in t he literature (See Table??). We fit the photometry to spectral tem- 
plates in the lPickled (119981 ) library, providing us with a robust and empirical approach to mea- 
suring the b olometric fiux -F r ot; (a verage errors of ~ 1.3%). General details of this method 



are gi ven in Ivan Belle et al.l ( 120081 ). although contrary to the approach in Ivan Belle et al. 



(120081 ) where interstellar reddening is a free variable in the SED fitting, this work sets the 
value of the interstellar extinction identical to zero since all our targets are nearby. This en- 
sures consistency across all targets in our survey. However, we note that this approach does 
not attempt to take into account any potential circumstellar reddening that may be present 
in individual systems. Thus, it is possible that previously published stellar parameters for 
stars may be slight ly different (within la ), even though they are based on largely the same 
data (e.g, GJ 581; Ivon Braun et al.ll201ll ). We were not able to measure a bolometric fiux 
for visual binary GJ 702A and GJ 702B because the photometry for these two objects is 
blended due to their close separation. Thu s for these two stars, we use the luminosity values 
computed from lEggenberger et al.l (l2008al ) using bolometric corrections. 
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With knowledge of 6'ld, and -Fbol, we calculate stellar effective temperature Teff via 
the rewritten Stefan-Boltzmann Law 



rEFF(K)= 2341(FBOL/^^D)^ (2) 

where -Frol is in units of 10"*^ erg cm~^ s~^ and 6'ld is in units of mas. The total luminosity 
given by L = Atk^Fbol is also derived for the stars in the sample. We present these results 
and associated errors in Table O 



3.2. Photometry, Metallicities, Masses, and X-ray Brightnesses 



One of the main contributions of this paper is to provide empirical relations linking 
observables, namely broad-band photometry and metallicity, to the measured physical prop- 
erties of late-type stars. Due to the fact that these objects are ver y bright in the infr ared, 
we are unable to get reliable photometry from the 2MASS Catalog (ICutri et al.ll2003al ). Al- 
ternatively, we list in Table [5] the collection of broad-band Johnson BVRIJHK photometry 
for all the program stars, along with the reference for each measurement. 

Metallicities for M-dwarfs are difficult to calibrate accurately due to their complex 
spectral features and molecular bands contaminating the continuum. [Fe/H] can be used 
for diagnostic purposes to understand the effects of a star's metallicity on relations between 
other stellar parameters. This is a crucial point of interest when comparing the stars to 
models, in which opacities in the atmospheres of M-dwarfs may not have properly been 
accounted for. The metallicity values for the stars in this paper were preferentially collected 
from a uniform source - so that we could minimize any systematic offsets with respect to 
different methods/references in the metallicity scales. Unfortunately, a single catalog of 
stellar metallicities does not contain all of the stars in this survey. 

For this study , we c ite the metallicities of M-dwarfs from the calibration presented in 



Roias-Avala et a. 



Neves et al 



(I20121) . wherever available. Alternatively, we quote the metallicity values 
J201lh. e.g.. for GJ 15 A and GJ 880. We find that all of the K-dwarfs 



from 

have entries in lAnderson &: Francis! (120111 ) , who provide average naetall icity values of several 
available literature references. While the lAnderson fc FrancisI (120111 ) metallicities are an 
average - and thus not of unif orm origin - we cross ch eck these averaged metallicities with 
those in the SPOCS catalog (jValenti fc Fisched l2005l ). where entries are available for the 
majority of our stars, and we find that the values are within 0.04—0.05 dex each other, 
revealing no systematic offset. The [Fe/H] values used in this study and references can be 
found in Table O 
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in order to compute stellar masses for the stars in the s ample, we emp l oy an absolute 
i^'-band mass-luminosity relation. While the relation from iDelfosse et al.l (120001 ) is com- 
monly used for M-type stars, it do es not extend to be valid f or K-type stars. Due to this, 
we chose to use the relation from iHenry fc McCarthy! (119931 ) because it is good through- 
ou t the whole rn a ss ran ge of the sample presented in this paper. We note that figure 1 



m 



tion derived in 



Delfosse et al.l (120001) illustrates c onsist ency in masses derived via the piece-wise solu- 



Henry fc McCarthyl (119931 ) and the polynomial fo rm of their solution in 



Delfosse et al. (2000). We apply the conversion from LeggettI ( 19921) to transform the John- 
son K magnitudes i n Tab le |5] into the CIT system (lElias et al.l Il982l ) so the relations in 
Henry &: McCarthyl (|l993l ) could be applied. These corrections are on the order of 0.02 mag- 
nitudes. The resulting estimates of the stellar mass are included in Table |6] and assume a 
10% uncertainty. 



Lastly, we use equation 1 in iLopez-Moraled (120071) to convert h ardness ratios from the 
ROSAT All-Sky Survey Bright Source Catalogue (IVoges et al.l 119991 ) to X-ray flux for each 
object and compute the X-ray to bolometric luminosity ratio Lx/Lbol, a useful diagnostic 
to the activity levels in late-type stars. These values are presented in the last column of 
Table El 



3.3. Literature Angular Diameters and Fundamental Properties 



At present, there are 10 published papers that present interferometricall y determined 
radii o f K and M dwarf star s with precision t o better than 5%: Lane et al. (I2OO1I): Segransan et al. 
(|2003 [ ): lBerger et al.l (l2006[) : Idi Folco et al.l (l2007h:lKervella fc FouauelJ2008h:lBovaiian et al. 



(I2OO8D : Ivan Belle fc von BraunI ( l2009f ): iDemorv et al.l ( l2009f ): Ivon Braun et al.l ( 1201 ll . 120121 ). 
These results within these papers amount to 22 measurements of 17 unique sources, and are 
compiled in Table [6l which also contains the values and references for the stellar metallicities 
(§133). 

We perform a SED fit for each of these stars to determine the bolometric flux, luminosity, 
and effective temperature in the same manner as we did for our targets in Section 13.11 The 
stellar mass is again derived from the /T-band mass-luminosity relation in lHenry fc McCarthy 
( I1993I ). Lastly, we also list the derived values of the derived bolometric to X-ray luminosity 
ratio Lx/Lbol- 
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3.4. Stars with Multiple Measurements 

Table [6] lists each individual measurement in our survey, as well as in the literature in the 
top two partitions. We also mark in Table [6] stars that have multiple measurements with a ^ or 
a The ^ marked stars include the stars discussed in Section [2] that we used for consistency 
checks across different instruments and interferometers: GJ 15A, GJ 166A, GJ 380, GJ 411, 
and GJ 699. Note that both GJ 380 and GJ 411 have two measurements in the literature. 
The marked stars indicate stars that have multiple literature measurements, but are not 
included in our survey: GJ 820A, GJ 820B, and GJ 887. We determine the weighted mean 
of these stars' radii and temperatures, and list them in the third portion of Table El The 
other properties of these stars, i.e., the spectral type, metallicity, bolometric flux, luminosity, 
mass, and Lx/Lbol, are not affected by the combination of multiple measurements, and we 
do not reprint them in this section of Table [6] for clarity. 

In summary, there are 17 late-type dwarfs in the literature with high-precision (better 
than 5%), direct measurements their radii, of which 10 are K-dwarfs and 7 are M-dwarfs. 
This publication adds equivalent (or higher precision) measurements of 16 additional late 
type stars (7 K-dwarfs and 9 M-dwarfs), plus 5 additional measurements of the 17 dwarfs 
in the literature for the purpose of consistency checks across facilities and techniques. The 
total number of directly measured radii of late-type dwarfs with precision better than 5% 
thus stands at 33: 17 K-dwarfs and 16 M-dwarfs. 



4. Color - Metallicity 



[Temperature, Radius, Luminosity] Relations for K- 
and M-Dwarfs 



In this Section, we determine relations between the observable properties of stars: broad- 
band color index and metallicity, to our empirically measured quantities: effective tempera- 
ture, radius, luminosity. To achieve this, the observables are fit to a multi-parameter, 2nd 
order polynomia l function, a.nd so lutions are found using the non-linear least-squares fitting 
routine MPFIT jMarkwardtll2009h . 



The following analysis uses our results plus the compilati on of results from th e literature, 
i.e., data presented in TableO We note that GJ 551 (M5.5 V. lHawley et al.lll996l ) is the latest 
spectral type star in the sample, departing from the rest of the bunch by ~ 1.7 magnitudes 
in its {V — K) color. This places it toward the extreme boundaries in the relations, and 
thus introduces a bias in the fit. Although we include it in our color- metallicity-temperature 
fits (Equation [3]), we caution that the properties of such late-type stars are not properly 
accounted for in this work. The coolest spectral type covered in these relations is M4 V. 
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As for the color-metallicity-radius, color-metallicity-luminosity relations (Equations H] and 
|5]), we do not include GJ 551 in the analysis because the adopted functional form of the 
relation is not appropriate over such a long baseline. This is covered in more detail in the 
sections to follow. The metallicity of the sample ranges from [Fe/H] = —0.68 to +0.35, 
where the median metallicity value is sub-solar at [Fe/H] = —0.19. The full sample of data 
is used in the fitting procedures regardless of metallicity, however we impose limits to where 
the relations hold true to where there is uniform sampling within the data set: stars with 
metallicity values ranging from —0.5 to ~ +0.1 dex. 

The specific form of each function (i.e., color- metallicity-temperature, color-metallicity- 
radius, color-metallicity-luminosity.) is defined in the sections to follow. For the temperature, 
radius, and luminosity relations, we find that adding the metallicity as an additional param- 
eter improved the fit up to a factor of two, where the strongest influence is exhibited on the 
redder (M-type) stars in the sample. Our choice of polynomial functions is based merely on 
the way the data appeared to present themselves, and has no physical relevance. For the 
radius and luminosity relations, the use of higher order (> 2) polynomials to fit the data 
was tested, but did not model the data any better. We stress that the functional form of 
the polynomials used for the radius and luminosity relations are strictly valid only for the 
range of stars specified. Beyond this range on both the hot and cool ends, the structure 
of the curve modelling the stellar properties will exhibit inflection points (see discussion in 
the sections to follow), requiring the use of higher order functions. We find that the major 
weakness in using a higher order polynomial is the lack of a more uniform sampling the data 
at the endpoints of the relations with respect to an object's color and metallicity, leading 
to poor constraints on the fits. For this reason, we do not incorporate these solutions in 
this work. However, we are able to ameliorate this aspect on the blue end of the fits by 
adding data for G-type stars to establish realistic (and empirical) boundary conditions in 
the color-metallicity-temperature relations, while still using a 2nd order polynomial. This is 
discussed more in Section 14.11 

The form of the function relating stellar temperature to other parameters varies in the 
literature. For instance, color-temperature relations (without metallicity) have been ex- 
pressed ranging from a 2nd order polynomial (e.g, see Blackwell fc Petfordlll99l[ ) to a 6th 



order polynomial (e.g., see iGrayl Il992l : iBoyajian et al.l |2012[ ) , and even a power function 



( Ivan Belle fc von BraunI 120091 ). When including the color and metallicity as variables, the 
preferred form is the one we adopt here, although it is parametrized t o contain the recip 



rocal temperature, 6'eff = 5040/Teff, as the fitting constant (e.g., see lAlonso et al.lll996l ). 
Additional corrections have been applied by fitting a 6th orde r polynomial to the residual s 
of the 2- variable color-metallicity function, e.g., in the work of iRamfrez fc MelendezI ( 120051 ). 
We anticipate that in the near future, new data published on the fundamental properties 
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of these types of stars will remedy this shortfalling and allow for robust testing to the true 
form of these relations (See Section [7]). We conclude this discussion by introducing a look-up 
table presenting the average temperature and radius for each spectral type in the sample 
(Section |0]). 



In each of the following sections defining the empirical color-metallicity-temperature 
(Section l4.ip , color- metalhcity-radius (Section l4.2p . and color- metallicity-luminosity (Sec- 
tion I4.3p relations, we overlay model predictions to compare with our results. However, we 
caution that a straightforward comparison of our solutions with models is indirect when 
regarded on the observational (color) plane, for the physical solutions of models' predictions 
must be transformed to the observational plane through color tables. Errors introduced 
i n color table conversions, whether synthet ic or empirical based, are still under scrutiny 



(iDotter et al.ll2008l : IVandenBerg et al.ll2010l ). Due to this, we refrain on discussing in great 
detail any incongruity within the observations and models on the color planes. In Section [5] 
however, we bypass the color planes and relate the physical quantities of stars directly. In 
these such cases, we are able to fully describe and quantify any anomalies in the model 
predictions compared to our observations. 



4.1. Color-Metallicity-Temperature 

In this Section, we derive empirical color-metallicity-temperature relations with the 
combined data set (Table [6]), using the indices of the Johnson photometric system: {B — V), 
{V - R), {V - I), {V - J), {V - H), and {V - K) (Table [5]). In order to properly model the 
relations at t he blue end of our d ata range, we add an additional seven G-type stars from 



our Paper I (IBoyajian et al.ll2012l ). whose properties were determined in the same manner 
as this work. We provide a table of the G-stars from Paper I in Table [71 While other G-type 
star angular diameters are available in Paper I, we find that only the seven listed in Table [7] 
have a complete photometric collection of all BVRIJHK magnitudes. So in order to ensure 
consistency of the improvement that the inclusion of additional measurements of G-stars will 
have on the color-metallicity-temperature relations we derive, we selected only these seven 
G-stars that will contribute equally in the relations for all color indices applied in this work. 

The solution for temperature as a function of color and metallicity is in the form: 



Teff(K) = ao + aiX + a2X'^ + a^iXY + a^Y + a^Y^ (3) 



where variable X is the color index and Y is the stellar metallicity [Fe/H]. The coefficients 
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ao-as and a statistical overview containing the number of data points used in the fit {n), the 
approximate color range over which the fit is valid, and the median absolute deviation about 
the best fit for each of the relations are presented in Table . The data and solutions are 
plotted in Figure Ei We color code the data plotted in these figures to reflect the metallicity 
of the star: blue indicating metal poor and red metal rich (a [Fe/H] range of —0.68 to +0.35). 
We also show the solution to the polynomial fit as colored solid lines illustrating a range of 
solutions for fixed metallicity values. 

All solutions give a median absolute deviation in temperature between 43 K and 70 K. 
The tightest of these relations are the \^-band to infrared color-indices {{V — H), and {V — 
K)), perhaps because they provide the most sensitivity as to where the peak of the stellar 
energy distribution is located, and hence a better determination of the effective temperature. 
However, we find that the {V — J) color gives the largest scatter. We are uncertain about 
the source of the deficiency in this bandpass index, but speculate that the systematics seen 
for 1.5 < {V — J) < 2.0, where all the points fall beneath the curve, and 1.2 < {V — J) < 1.5, 
where the points all fall above the curve, is indicitave of the fact that a higher order function 
might be needed to properly model the rapidly changing slope of the data. 

In Figure [6] we show our solution as iso-metallicity lines for fixed values of metallicity 
in 0.25 dex increments. It can be seen that each of these iso-metallicity curves converge for 
the earlier-type stars. We find that adding metallicity as an extra parameter in the color- 
metallicity-temperature relations only improves the median absolute deviation of the data 
to the fit by ~ 30%, with the most prominent effect between 3500 and 4000 K, where the 
solution for [Fe/H] = —0.5 is about 100 K cooler than a solar abundance of [Fe/H] = 0. 
This is true even for the color-indices with long baselines, (i.e., {V — J), (V — H), and 
{V — K)), which we find to be just as sensitive to metallicity as the shorter baseline color- 
indices examined here. This could be caused by the redder apex of the SED in these types 
of stars and the formation of molecular features in the atmosphere that contribute to larger 
flux variations in the IR band fluxes. Although this is a weak detection to the influence of 
the metallicity to color and temperature, it can perhaps be refined with a better sampled 
range of metallicities, especially for high and low values of metallicity at the cool endpoint 
of the relations. 

We include the {B — V) color-metallicity-temperature relation but caution against using 
this relation for stars with peculiar abundances since {B — V) color is known to be strongly 
affected by stellar metallicity. In fact. Figure [6] shows one of the most extreme outliers 



''Note that some of the K and M stars do not have pubUshed measurements in aU photometric bands, 
and can be seen when the total number of points used in the fit n is less than 33. 
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toward the bluest of {B — V) colors and several sigma under the curve is the metal-poor star 
GJ 53 A (yU Gas A). We note that GJ 53A does not fall on the (5 — \^)-metallicity-temperature 
relation, but it is not an outlier in any of the other color-metallicity-temperature relations 
(see Figure E]). We chose to include it in our color-metallicity-temperature analysis to be 
complete, and we note that it does not change our results when removing it and re-fitting 
the data. In fact, all solutions except for the {B — V) color-metallicity-temperature relation 
are shown to converge at temperatures hotter than ~ 4000 K. 

In Table [S] we show that the median absolute deviation of the fits using the extended 
GKM sample is within 5 K of the KM sample. This is shown graphically in Figure O where 
we plot the G-dwarfs used to constrain the fit on the hotter side, along with the K- and 
M-dwarf sample, showing that there is a smooth transition of temperature extending to the 
bluer colors. Thus, the extrapolation of the solutions on the hotter side is valid through 
spectral type GO. 

On the other hand, extrapolation to stars redder than the ranges specified in Table [8] 
is not recommended. Figure shows that past this specified range, there is a large gap in 
color index before the reddest endpoint, the star GJ 551, is encountered. So, although we 
use the data for GJ 551 in the fitting procedure, the paucity of data in this range makes 
the results unreliable for such red colors, as for the relations presented here, the position of 
GJ 551 deviates several hundred Kelvin from the modelled curve. 

In Figure E] we show the — fC)-temperature solut ion for dwarfs (expressed as a p ower 
function) derived from interferometric measurements in Ivan Belle fc von BraunI (120091 ). dis- 
playing excellent agreement with our fit. Also for comparison, we show the model based 
temperature scale for the s olar metallicity (-B -- V^), (V — J), (V — H), and (V — K) color- 
temperature solution from iLejeune et al.l ( 1l998l ) in Figure O While the temperatures from 
Lejeune et al.l ( 1l998l ) for the earlier type stars show excellent agreement with our own, the 
{y — J), {V — H), and {V — K) solution diverges at temperatures < 4500 K (approxi- 
mately the K5 spectral type) to produce temperatures ~ 200 K elvin cooler than the ones 
that we measure with interferometry. For the {B — V) colors, the iLejeune et al. Jl998h scale 
diverges to predict higher temperatures from 5000 < Trff < 4000 K . Belo w this tempera- 
ture, the data are modelled more satisfactorily with the ILejeune et al.l (119981 ) scale, although 
with a slope much different than our own. The difference in slope is no surprise, for the 
low-order polynomial we use in this work is not capable of modelling any rapidly changing 
slope and/or kinks in the data. Lastly, we sh ow the BT-Settl PHOENIX model atmosphere 
color-temperature curves (lAUard et al.ll2012l ) for solar metallicity in the (5 — V"), {V — J), 
and {y — K) plots, expressing by far the best agreement with the temperatures we derive 
here for all ranges of temperatures. It is only in the {B — y)-temperature solution that for 
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stellar effective temperatures less than 4000 K, the BT-Settl model predicts a steep drop in 
temperature at colors too red by about [B — V) ~ 0.2 magnitudes. This is due to troubles in 
de termining the B rn agnitudes from the synthetic PHOENIX model grid (see the discussion 



m 



Potter et al.l 120081 ). 



Specific references citing temperatures of these stars in Table [6] are plentiful. For in- 
stance, for any given target in the sample, there are, on average, tens of references to the 
effective temperature in the literature, estimated using a number of different techniques. As 
such, a quantitative comparison of our temperature values to all temperature measurements 
in the literature is beyond the scope of this paper. Instead, for each star we prepare an 
overall qualitative assessment of how our temperatures compare to previously obtained ones 
by examining the full range of temperature estimates listed in Viziei^l. Figure [7] shows the 
results of this exercise, where the x-errors indicate our measurement precision, and the y- 
errors indicate the range in temperatures found for each object. Note that a search in Vizier 
may not gather all references to the temperature in the literature, but rather a collection of 
available on-line material. In fact, a search in Vizier resulted in no references for tempera- 
ture estimates of GJ 702B, and only one for GJ 551. We thus performed a more thorough 
literature search to gather a range of temp erature estimates for thes e sta rs to make this com- 
parison (the addi tional references include Eggenberger et al. 2008a and Luck fc Heiter 20051 



for GJ 702B, and Morales et al.ll2008l : IPemorv et al 



2009 



Wright et al 



12011 



for GJ 551). 



4.2. Color-Metallicity-Radius 

Similar to §4.H we use the data in Table to derive empirical color-metallicity-radius 
relations in the form of Equation HJ where X is the color index, and Y is the stellar metallicity 
[Fe/H]. 



R{Rq) = ao + aiX + a2X^ + asXY + a^Y + a^Y^ (4) 

However, as opposed to the temperature relations, we refrain from including the G-type 
type stars in Paper I to constrain the bluer end of the relations, as their radius may be 
influenced by stellar evolution (we assume that a star with radius less than 0.8 has not 
had time to evolve off the main-sequence). The inclusion of GJ 551 is also omitted from this 



^Vizier provides access to the most complete librar y of published astr o nomic al catalogues and data tables 
available on-line. For a description of the catalog, see Ochsenbein et al. ( 2000f) 
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analysis, for the adopted form of the equation we use does not allow for extrapolation of the 
curve to such red colors. 

The coefficients and statistical overview of the solutions are in Table M We find that, 
when adding metallicity as an extra parameter adds significant improvement to the signif- 
icance of our fits, cutting the median absolute deviation of the observed versus calculated 
radius in half. Figure [H] displays the data and our fit using the same color scheme as identified 
in the temperature relations to illustrate the metallicity. The data show a large spread in 
radii for a given color index toward the coolest stars (e.g., for (B — V) ~ 1.5 and (V — K) ~ 4, 
the radii range from 0.4 — O.6-R0). We discuss this effect in more detail in Section [531 

The ranges of color indices where the relations are valid are listed in Table [91 The ranges 
listed are hard limits, and do not allow for extrapolation. There are several reasons for these 
strict boundaries. On the blue end, we are not able to observationally constrain the form of 
the upper part of this curve due to stellar evolution. This causes the unphysical behavior of 
the curve to 1) fold over as a parabola, and 2) display cross-overs among fiducials of different 
metallicities (see Figure[H]). In actuality, a zero-age main sequence radius curve would have an 
inflection point near this upper boundary we impose on the ranges in Table [HI Past this point, 
the stellar radius will continue to rise with decreasing color index (but unfortunately also 
changes with increasing age). The use of the color-metallicity-radius relation (Equation H]) 
on the blue end, despite the cross-over, will still lead to results within the stated errors, so 
long as the object in question is non-evolved. The red ends of the relations are strictly hard 
limits as well, where, by visual inspection of Figure [HI extrapolation would result in null radii 
estimates approaching the late M-type star regime. 



synthetic color transformations (lBesselllll990l : ICutri et al. 



We show 5 Gyr Dartmouth i sochrones (iDotter et al.ll2008ll ba sed on the U BV{RI)cJHKs 



2003al ) within the panels display- 



ing the {B — V), (V — J), {y — H), and {V — K) data and relations in Figure [HEI- These 
isochrones are for two metallicities, [Fe/H]= and —0.5, roughly encompassing the metal- 
licity range of the data included in this paper. The {B — y)-radius plot also shows the 
Dartmouth isochrones under the same conditioii s, but with the semi-empirical BV{RI)c 



color-transformations (IVandenBerg &: ClemI [20031 . blue lines). Dartmouth models for John- 



son R and / colors are not provided, so we are unable to compare these color indices to our 
observations. 

There is acceptable agreement of the models and data with the earlier type stars, and 
the isochrones also show the sharp change in slope to increasing radii when extending beyond 



^Since conversions between K and Kg of these stars are on the order of the magnitude errors, we assume 
Ks ~ K in these plots. 
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our data set to bluer colors (G-type stars). Extending to the later-type stars however, we 
find that for stars with radii around O.6R0 and below, the models predict a drop in radius 
at colors too blue by several tenths of a magnitude compared to what we observe in the 
{V — J), {V — H), and {V — K) colors. The {B — V) panel shows that at this point, 
the solar metallicity model unexpectedly drops and crosses over the model —0.5 dex curve. 
This is a well documented failure of the synthetic c olor grid's ability t o reproduce accurate 



colors for short wavelengths (e.g., see discussion in iDotter et al.ll2008l for the accuracies in 
the color conversions). For instance, the use of the semi-empirical color transformations in 
the {B — V) panel (blue lines) do not show this fiaw, demonstrating the contrast of using 
different color-table conversions. 

Overall, the models show a pattern of divergence with the observations in the color- 
radius plane. The reason behind this disagreement is not easily identified, for it could be a 
consequence of under-predicted radii for a given color index, or an offset in the color-tables 
used to calculate the colors of the later- type stars (see Section [5]). 



4.3. Color-Metallicity-Luminosity 

In this Section we use the data in Table O to derive relations between a star's color and 
metallicity to its luminosity, which can be inverted to solve for the distance, if unknown, so 
long as one can assume that the star is on the main sequence. In the same way as the radius 
and temperature relations, we express luminosity as a function of color and metallicity using 
a polynomial in the form of: 



logL(L0) = ao + aiX + a2X^ + a3XY + a^Y + a^Y^ , (5) 

where the variable X is the color index and Y is the stellar metallicity [Fe/H]. Table [TO] lists 
the coefficients and overview of each relation, and Figure M show the data and solutions for 
our fits. Similar to the radius relations in Section we find that the addition of metallicity 
as an extra parameter reduces the median absolute deviation by a factor of two. The plots 
in Figure [9] show that the solutions converge for the earlier type stars regardless of the 
metallicity. However, for the later-type stars there is a spread in luminosity for a given color 
index, confirming that the metallicity is a necessary factor in estimating a star's luminosity 
using color relations. 

We impose the same hard limits for the range of color indices over which the color- 
metallicity-luminosity relations are valid, just as we do in the radius relations. That is, no 
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extrapolation of the curve is warranted for stars bluer or redder than the indicated color 
ranges listed in Table [TOl Once again there is a tendency for some of the iso-metallicity 
solutions displayed in Figure M to cross over each other. This effect is non-physical, and is 
a product of the metallicity dependent solutions converging at the blue end of the relations. 
Despite the appearance of the cross over, any application of the relations with real data will 
still lead to reliable results with the quoted errors - so long as the boundary conditions are 
carefully regarded (Table [T0|) . 

The luminosity versus color plots can be interpreted as empirical isochrones. Thus, the 
panels displaying the ( B — V), (V — J), (V — H), and {V — K) relations in Figure [9] also show 



Dartmouth isochrones (iDotter et al.ll2008l ) for an age of 5 Gyr and two metallicities 0.0 and 
—0.5, roughly bracketing the metallicity range of the data included in this paper. Within the 
{B — V) panel, the Dartmouth model solutions using both the synthetic (black) and semi- 
empirical (blue) color conversions are shown, as described above in Section 14.21 Two items 
are immediately clear based on the comparison of the models to our data. First, the models 
reproduce the properties of the bluer stars (i.e., stars with logL/L© > —1). Below this limit, 
the models predict an observed steep drop in luminosity, however, this drop is predicted at 
colors bluer than what we observe. This is very similar to our color-radius results (Section l4.2( 
Table [9l Figure [8]), and provides substantial motivation for modellers to identify whether or 
not the source of this disagreement is within the computed astrophysical properties, or in 
the tabulated colors of the later-type stars. We discuss this more in Section |5l 



4.4. Relations to Spectral Type 

For each spectral type observed, we use the data in Table in order to build a look-up 
table referring a spectral type to the star's temperature and radius. These results are in 
Table [TTl The quantity n indicates the number of stars that fall into that spectral type bin. 
The value of the parameter given is the average value of all stars within the spectral type 
bin, and the a is the standard deviation of the parameter uncertainties for each spectral type 
bin. The spectral types with only one measurement (n = 1) simply lists the individual value 
and the measured error of that measurement. 
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5. Relations to the Global Properties [Luminosity, Temperature, Radius, 

Mass] for K- and M-Dwarfs 

The Sections I5.H I5.2[ I5.3[ and 15.41 link the global stellar properties to each other: 
luminosity-temperature, luminosity-radius, temperature-radius, and mass-radius. In each 
Section, we derive empirical relations based on the data of the 33 stars in Table El Following 
the observational results, we introduce and discuss model predictions of each star, as well as 
compare overall agreement to our relations and model relations of the physical quantities. 

In the following analysis, we use the Dartmouth models in the interpretation of our 
data because of the flexibility and user-friendly interface of the web-based interpolator that 
generates customized isochrones for user^. We note that there are several evolutionary 
models for low-mass stars available to the community. F igure [TO] shows our d ata and four 



models: Padova, Dartmouth, BCAH, and Yo nsei-Yale (iGirardi et al.l l2000t iDotter et al. 



20081 : iBaraffe et al.lllQQSt iDemarque et al.ll2004l ). In the luminosity-temperature, luminosity- 



radius, and temperature-radius plots, the dashed line corresponds to the solar metallicity 
model (see color legend within plots), and the hashed region spans a metallicity of to 
—0.5 dex. In the mass-radius plot, the solid and dashed lines mark a change in metallicity 
from [Fe/H] = to —0.5. We restrict the metallicity range plotted to conserve consistency 
across models (BCAH98 models only cover this range in metallicity). Note that for the 
BCAH98 models we use mixing length of a = 1.0 for stars < 0.6 Mq and the solar calibrated 
mixing length of a = 1.9 for stars > 0.6 M0. Also, only a solar metallicity solution with 
the BCAH98 models is shown for the stars > 0.6 Mq for it is the only model available. Of 
these four models, the Dartmouth and BCAH show to reproduce the trends of our data the 
best. While each model is approximately consistent for stars brighter than 0.16 Lq (earlier 
than ~K5), the conclusions here are speciflc to the Dartmouth models, and should reflect 
the best case scenario to the observed offset with models for stars below this mark where 
the agreement is poor across model predictions (See Figure [TO|l . 



5.1. Stellar Luminosity versus Temperature 

We use the 33 stars in Table [H] to derive a relation between the luminosity as a function 
of temperature expressed as: 



^http://stcllar.dartmouth.edu/niodels/isolf.html 



- 21 - 



logL(L0) 



5960.5710(±207.7541) + 4831.6912(±172.6408)X 
1306.9966(±47.8104)X^ + 117.9716(±4.4125)X2 



(6) 



where the variable X stands for the logarithmic effective temperature log(TEFF)5 and is valid 
for non-evolved objects within a temperature range of 3200 to 5500 K. A metallicity depen- 
dent solution did not show any improvement to the fit, therefore this relation is metallicity 
independent. The median absolute deviation of this fit is 0.0057 L©, similar to the abun- 
dance dependent color- metallicity-luminosity relations (Equation [5l Table [TO|l . The top plot 
in Figure [TT] shows the data and the solution to this fit. The middle plot in Figure [TT] 
shows the fractional deviation of observed luminosities compared to the empirical relation 
expressed in Equation [61 

The fractional residuals in the middle panel of Figure [TT] (calculated as (Lobs— -^Fit) / L-pn, 
where Lp^ is the luminosity of the empirical relation at the star's observed temperature) show 
that for stars with temperatures above ~ 3500 K, the observed luminosities lie within ±10% 
of our empirical solution. However, we note that for stars with temperatures below ~ 3500 K, 
the scatter about the fit is much higher. The plotted residuals also show no pattern with 
respect to the stellar metallicity. This reinforces the fact that we were not able to solve for 
a metallicity dependent solution for Equation [6] using our data. 

In the top plot, we overlay Dartmouth 5 Gyr isochrones for [Fe/H]= 0, —0.5 (dash- 
dotted and dotted lines, respectively). We find that the observed agreement of our data 
and the Dartmouth models in Figure [TT] is fairly decent, although much more favorable in 
the earlier type stars. To evaluate this agreement of our data and the Dartmouth models 
more quantitatively, we generated Dartmouth model isochrones for each of our stars at its 
respective metallicity. Each isochrone was chosen at an age of 5 Gyr (appropriate for non- 
evolved low-mass field stars) with the default Helium mass fraction {Y = 0.245 + 1.5 * Z), 
and alpha-enhanced elements scaled to solar ([a/Fe]= 0)§. We then interpolate through 
the model grid to find the model temperature, TMod, at the observed luminosity of each 
star. The results are plotted in the bottom panel of Figure [TT] The results show that for a 
given luminosity, models overestimate the stellar temperature by a average of 3% (~ 120 K) 
for stars with temperatures under 5000 K. For the 7 stars with temperatures greater than 
5000 K, we detect a trend that the temperatures of the more metal-rich stars (orange points; 
[Fe/H] ~ 0.0) are underestimated by the models by a few percent. On the other hand, 
the more metal-poor stars (green points; [Fe/H] ~ —0.25) with temperatures greater than 



'We use a [a/Fe]= 0.2 for GJ 53A (See discussion in 



Bovaiian et al 



20081) . 
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5000 K show that their temperatures are overestimated by the models by a few percent. 
This effect is also seen for stars cooler than 5000 K, where the less abundant the star is in 
iron, the more offset the star's predicted temperature will be with the observed temperature. 
This pattern of metallicity within the residuals in the bottom panel of Figure [TT] is telling 
that the models could be over-predicting the influence of metallicity on the global stellar 
parameters. 



We use the 33 stars in Table [HI to derive a solve for an solution on the radius-luminosity 
plane expressed as: 



where X is the stellar radius in solar units. The median absolute deviation about the fit is 
0.008 Lq, and the data and solution are plotted in the top panel of Figure [T2l The middle 
panel in Figure fT2l shows the fractional deviation in luminosity about the fit. Stars with radii 
between ~ 0.4 Rq and ~ 0.67 show the tightest correlation, deviating from the fit by 
less than ±10%. Outside this range in radii, the scatter is worse, where on the lower end the 
scatter increases to double this level and on the upper end the scatter increases to triple. We 
observe no pattern with metallicity in the fractional residuals compared to our fit, enforcing 
the fact that a metallicity dependent solution is not appropriate on the luminosity - radius 
plane. 

We plot Dartmouth 5 Gyr isochrones for [Fe/H]= 0, —0.5 on the top portion of Figure [T2l 
By inspection of these curves, the solar metallicity model agrees well with our empirical fit 
(Equation [TJ while the metal-poor isochrone appears slightly offset from the data, similar to 
the luminosity - temperature plots in Figure [TT] To make a fair assessment of the comparison 
to each star's unique composition, we use the custom-tailored Dartmouth isochrones gener- 
ated for the 33 stars (discussed above in Section [5TT]) and interpolate through the model grid 
to determine the model radius -RmocI at the observed luminosity for each star. The fractional 
deviations of the observed radii from the model radii are presented in the bottom panel of 
Figure [T2J This shows a clear offset on the order of 5% for models to under-predict the 
stellar radii at a given observed luminosity. 



5.2. Stellar Luminosity versus Radius 



logL(L0) 



3.5822(±0.0219) + 6.8639(±0.1562)X 
7.1850(±0.3398)X2 + 4.5169(±0.2265)X3 



(7) 



The deviations in observed radii to those predicted by models become prevalent for 
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stars with radii < 0.7 R© (bottom panel of Figure [T2l) . For the full range of stars in our 
sample however, we also detect the same pattern in the residuals shown in the bottom panel 
of Figure [TT|, where the stars with lower metallicities have radii that are are more deviant 
than model predictions than the stars with more solar type metallicity values. While the 
effect is less pronounced on the radius plane than on the temperature plane, the conclusions 
remain the same: models appear to be over-predicting the impact stellar metallicity has on 
the physical stellar properties. 



Our results are consistent with the findings in the literature (e.g. |Poppeiill997l : lRibas et al 



20071 : lL6pez-Moralesll2007l : [Morales et al.ll2008[ ). where models are shown to predict tempera- 
tures too high and radii too small while still being able to correctly reproduce the luminosity. 
Such references however, are typically referring to binary star properties, and we note that 
this result is quite novel to detect such a distinct pattern for single stars. Additionally, a 
unique contribution to our project's analysis is that we are able to reveal the pattern of the 
offset in connection with the stellar metallicity, an observable often-times unreachable for 
binary stars due to the complexity of the blended spectrum from each component. 



Casagrande et al.l ( l2008l ) use the IRFM (infrared-flux method) in order to determine 
properties of M-dwarfs. Their conclusions show that there is a radii offset present when 
comparing their observations of single stars to models. However, while similar to our own 
conclusions, they estimate this offset is on the order of 15-20%, unlike the 5% we claim 
from our analysis. This difference is possibly caused by the errors introduced in their semi- 
empirical determinations of the stellar properties. They also identify a discontinuity at the 
transition from K- to M-dwarfs, most apparent on the luminosity - temperature plane, that 
is not reproduced by models. We do not observe this feature, however, this region from 
4200 to 4300 K is not well populated within our interferometrically observed sample, as we 
only have one data point that lies in this rang e of temperatures. Perh aps what is more 
interesting is a region not explicitly identified in ICasagrande et al.l ( 120081 ) ; the point where 
the deviations of their data compared to models begins. The data in their figure 12 show 
this is at ~ 5000 K (Mbol ~ 6.5). Their figure 13 (a) and (c) show this point at ~ 0.7 R© 
and ~ 0.7 Mq, the exact positions we identify the beginning of the trend in our data. 



5.3. Stellar Temperature versus Radius 

In order to link effective temperature to stellar radius, we use the 33 stars in Table [6] to 
derive a single parameter solution in the form of a 3rd order polynomial, valid for non-evolved 
objects with temperatures between 3200 and 5500 K: 
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R{Rq) = -10.8828(±0.1355) + 7.18727(±0.09468) x lO^^X 

-1.50957(±0.02155) x lO^^X^ + 1.07572(±0.01599) x 10"^°X^ (8) 

where the variable X is the effective temperature. Like the radius - color plane ( §4.2p . there 
is a lot of scatter in the radius of a star for a given temperature (see Figure [T3!) . However, the 
data on the temperature - radius plane does not show the same spread in the measurements 
due to differences in metallicity as in the color-metallicity-radius relations, as its solution 
is independent of metallicity. We find that the median absolute deviation of the observed 
versus calculated radius for Equation [8] is 0.031 R©. This value is similar to ones derived via 
the abundance dependent color-metallicity-radius relations of Equation H] shown in Table [91 

Stars more massive than the sample of K- and M-dwarfs (i.e., masses > 0.8 Mq) are 
sensitive to evolutionary effects. For example, with increasing age, a more massive star 
will first begin to evolve off its zero-age main-sequence position to larger radii and hotter 
temperatures. We illustrate this in Figure [T31 showing the expected model predictions of 
a 0.8 M0 star (blue *), a 0.9 M0 star (red +), and a 1.0 M0 star (black x) at the age of 
1 Gyr and 4.5 Gyr from solar-metallicity Dartmouth model isochrones. Arrows connecting 
these positions point to the top-left of the plot, i.e., the direction of evolution to larger radii 
and hotter temperatures. We also show in Figure [13] the Sun's present position. To model 
the data beyond our sample of K and M dwarfs to hotter temperatures, we re-fit the data 
adding the Sun as a point of reference to more massive stars, effectively making a solar-age 
calibrated isochrone. This relation is expressed as: 

R(Rq) = -8.133(±0.226) + 5.09342(±0. 16745) X 10-^X 

-9.86602(±0.40672) x lO^^X^ + 6.47963(±0.32429) x lO'^X^ (9) 

where X is the effective temperature. We show this solution in Figure [13] as a blue line, as 
well as the Sun's position (farthest left point in Figure [T3l assuming Tq = 5778 K, R = 1 Rq). 
The extrapolation of this curve past 5500 K is shown as a dotted hne that intersects the 
Sun's present position. A comparison of the temperature - radius relation excluding the Sun 
(Equation [8]) and including the Sun (Equation [9]) shows that while the hottest of the K-stars 
in our sample are modelled better by Equation [HI there is an inflection point to larger radii 
beyond our data sample range of spectral type KO to meet the Sun's position at its age 
today. 

What is more intriguing, is whether or not the models are able to reproduce the trend to 
the temperatures and radii of single stars, with the omission of the luminosity as a constant 
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like structured in the last two Sections of 15.11 and 15.21 The bottom-left panel in Figure [TO] 
shows that the Dartmouth models predict the sharp drop in stellar radii to happen at 
hotter temperatures than what we observe. This is also seen in the temperature-luminosity 
(Section 15. ip and radius-luminosity (Section 15. 2p relations displayed in Figure [TT] and [T^ 
where compared to our observations, the Dartmouth models predict temperatures too high 
and radii too low for a given luminosity. 

Similar to the method described in Sections 15.11 and 15. 2[ we use the custom-tailored 
Dartmouth isochrones for each star to evaluate expected parameters from the model com- 
pared to our observations. First, we assume that the temperature is a constant within the 
observations and model grid. For each star, we then interpolate within the model grid for 
a model radius at the observed temperature. The second approach is the reverse, and we 
assume that the radius is constant within the observations and model grid. A model tem- 
perature is then solved for by interpolation through the model grid at the observed radius. 
We show the results in Figure UM 

The top plot in Figure [H] shows the fractional deviation of temperature with the ob- 
servations to the Dartmouth models as a function of temperature. The bottom plot shows 
the fractional deviation of the observed to Dartmouth model radii as a function of radius. 
The fractional temperature offsets for stars with temperatures less than 5000 K, show that 
models predict temperatures higher than we observe by an average of 6%. The fractional 
radius offsets for stars with radii less than O.TR©, show that models predict radii smaller 
than we observe by 10%, where this offset increases with decreasing radii to ~ 50% for stars 
with radii ~ O.4R0. The three stars with temperatures below 3300 K (GJ 725B, GJ 551, 
and GJ 699) produce null results when interpolating at the observed temperatures to derive 
model radii, and thus their results are not plotted in the bottom plot of Figure [TH The com- 
puted offsets in temperature and radius given here reflect doubly deteriorating conditions 
compared to the results in Sections 15.11 and 15. 2[ where the stellar luminosities were held as 
a constant. These compounded errors should be regarded with caution when using models 
to relate stellar temperatures and radii for stars < 5000 K and < O.TR©. 



5.4. Stellar Mass versus Radius 



Data from double-lined spectroscopic echpsing binary (EB) systems are standards for the 
mass-radius relations available tod ay. For the fo ll owing discussion, we collect the binary star 
parameters presented in table 2 of [Torres et al.l (120101) togeth er with the low- mass binaries 
and secondaries in binaries compiled in iLopez-Moraled ( 120071 ). We impose the same criteria 
as for the single stars in our analysis, limiting the sample to only allow stars with mass 
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and radius < 0.9 Mq and R©, and where the radius is measured to better than a 5% error, 
leaving a total of 24 individual stars. These properties of the single stars are discussed in 
Section [31 where we have a total of 33 stars in this range of masses and radii. Note that 
although we limited the literature measurements of single stars to only include those with 
angular diameters measured to better than 5%, the sample remains unchanged if we filter 
with respect to the error in the linear radii to be better than 5%, due to the close proximity 
of the objects (within ~ 10 pc) and thus well known distances. 

To determine a mass-radius relation, we use a second order polynomial to fit the data 
for single stars (Equation [TOj) . We experimented with higher order polynomials and found 
that they did not did not improve the fit. We solve for a solution for the EB's in the same 
manner (Equation [TTj) . Data and fits are displayed in Figure [13 The resulting relations are: 



Rss(i?0) = 0.0906(±0.0027) + 0.6063(±0.0153)M, 

+0.3200(±0.0165)M^2 (10) 



Reb{Rq) = 0.0135(±0.0070) + 1.0718(±0.0373)M, 

-0.1297(±0.0367)M2 (11) 

where is the stellar mass in solar units. 

In Figure [151 we show the data and solutions for single (circles, solid line) and binary 
stars (squares, dotted line). These relations are both consistent with a 1:1 relation to the 
mass and radius of a star (dashed line). One can see from the figure that the radii of single 
and binary stars are indistinguishable for a given mass. This shows that tidal influences on 
binary components radii may not be a concern when viewed on this scale. In other words, the 
scatter in the measurements in EB systems may wash out the effect the binary period might 



have in the radii of the binary star components (for instance, see discussion in iKraus et al. 



201l|). 



5.5. The Endless Discussion About Model Predictions of Late- Type Stellar 

Masses and Radii Briefly Continues Here 



Historically, large discrepancies have been observed when comparing observed radii 
with radii predicted by models for low-mass stars (for example, see discussion in SJT] and 
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Lopez-Moralea 120071 ) in the sense that the models tend to underestimate the stellar radii 



at a given mass. We compare our measured radii to the model ra dii from the 5 Gyr 



model isochrones of Padova, Dartmouth, BCAH98, and Yq nsei-Yale (iGirardi et al.l 12000 



Potter et al.l l2008l : iBaraffe et al.l Il998l : IPemarque et al.l 12004 ) in Figure HDl Points for the 



masses and radii from eclipsing binaries are also plotted. 

Due to the density of information from showing several model predictions along with 
the data and errors, the visibility of the claimed radius discrepancy with models in Figure HU] 
is difficult to see clearly. While it is abundantly documented in the literature that such a 
discrepancy in the predicted and observed radii exists for binary stars, the intentions of this 
work is to show its equivalent - if present - for single stars. Already in Sections 15. II and 1 5.21 we 
have shown proof of this discrepancy in the luminosity - temperature and luminosity - radius 
plane, where models over-predict temperatures by ~ 3%, and under-predict radii by ~ 5% 
compared to our observations. These results presented on the luminosity - temperature and 
luminosity - radius plane are more robust for single stars because temperatures and radii 
are properties of single stars we can directly measure to high precision with interferometry. 
However, the ability to comprehensively show the radius discrepancy from this data set is 
quite difficult on the mass - radius plane, since the mass errors for single stars are character- 
istically large (as shown on the bottom axis of the plot). To first order however, we find it 
useful to present Figure [HI a series of plots showing the fractional deviation of the observed 
to theoretical radii for single stars. 

The theoretical radii, RmocI, are computed by using the mass (calculated from the mass- 
luminosity relation; Table |6]) and interpolating through the custom-tailored mass-radius grid 
of values computed for each star's metallicity in a 5 Gyr Dartmouth model isochrone. The 
result is a value for a model radius at the given mass and metallicity. The difference is 
then between the observed radius (from interferometry; Table E]) and the model radius (from 
interpolation of the mass from the M-L relation in the model mass-radius grid). Since the 
single star mass errors are quite large, the errors in the theoretical radii dominate the y-errors 
shown in Figure [161 Again, we caution that the absence of any signal in a (-Robserved — -Rmodci) 
plot could be misleading when using single stars as clues into this discrepancy on the mass- 
radius plane due to the large mass errors. 

On the top panel of Figure [161 "W'e show the difference in observed versus model radii 
as a function of mass. The two middle and two bottom panels show the fractional deviation 
of radius as a function of metallicity [Fe/H] and stellar activity level Lx/L^ol (Table [6l), 
since previously proposed diagnostics to explain any disagreement with models have been 
linked to abundances and stellar activity levels for stars in this mass regime. The bottom 
four panels are segregated to show stars of all masses < O.DM© and only admit stars with 
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mass < O.6M0 (left panel and right panel, respectively). 



The top panel on the plot shows that even with the stated large errors, stars with 
masses below 0.42 Mq (8 total) have bigger radii than those of the Dartmouth models, 
where 6 of the 8 deviate by several sigma. We identify these 8 stars with masses below 
0.42 M0 as blue points within the right hand side plots of Figure [T6l Regarding the fractional 
deviation in radii for the blue points with respect to metallicity, we find that for metallicities 
[Fe/H] < —0.35, 4 of the 5 blue points deviate by more than la. Any pattern with respect to 
Lx I -^BOL is less clear within the data. The fact that we do not see any trend with activity is 
not a surprise, due to the low levels of these single field stars compared to their active binary 
sta r counterparts (typic ally a decrease in X-ray flux by two to three orders of magnitude; 



see 



Lopez- Moralesll2007h 



Our conclusion supports the known discrepancy in observed radii compared to models 
for stars less than 0.42 M©, showing a slight dependence on metallicity. This conclusion is dif- 
ferent fro m previous works th at focused on interferometric radii of single stars. For instance, 
although iBerger et al.l ( l2006l ) show the presence of the radius discrepancy with models, they 
observe an opposite trend compared to our own: t hat stars with the m ost deviant stellar 
radii have larger metallicities. Contradictory t o the Berger et al.l (2006 ) conclusion as well 



as the one presented here is the analysis from iDemory et al 



(120091 ) ■ where they find that 



the radii of single stars are consistent with predictions of stellar models. 

The differences in the previously published results and the one presented h ere are multi- 



fold. A primary difference is the sample used in the analysis of each work. iBerger et al. 



(120061 ) applied no filter to the diameter precision allowed in their analysis, and at the time of 
publication, this amounted to 10 interferometric radii of M-dwarfs available for their analysis. 
Minding our filter of using only dat a with sizes determin ed to better than 5%, this constitutes 
to only 5 of the 10 stars used in the IBerger et al.l (120061 ). In addition to limiting the precision 
of measurements allowed in our analysis, we also tripled the sample size of M-dwarfs (14), 
making for better statistics and thus giving more weight to our results. Our work also uses 
recently defined metallicity calibrations for M-dwarfs (see references in Table [6]). These 
new references use techniques to measure metallicity that are mor e refin ed than the ones 
used to derive the metallicities for the stars cited in IBerger et al.l (120061 ) . and thus could 
explain the different conclu sions presented here. However, by comparing the metallicities 
used in IBerger et al.l (120061 ) to the ones we use in this work, we find an agreement within 
0.02 — 0.14 dex. This is less than the l-cr uncertainty in metallicity (~ 0.2 dex) and thus 
should not be the cause of new trend we observe. 



The sample of M-dwarfs at the time of the 



Demory et al 



(I2OO9I ) publication was identical 
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to the time of iBerger et al. J2006l i^. In their analysis, they include b oth K- and M-dwar f 
measurements, choosing to disregard all but one measurement of the iBerger et al.l (|2006[ ) 
data, and limit published measurements to have a precision better than 10%. 



The comparison of the radii offsets examined in iDemory et al.l (120091 ) come from an 
analysis of an expanded mass range encompassing both K- and M-dwarfs, similar to the one 
presented here on t he LHS panels of ou r Figure [161 In the full mass range, we come to the 
same conclusion as iDemory et al.l (120091) in their figure 9; that we see no deviation in single 
star radii with respect to metallicityl^°l. However, we find that by adding the stellar mass 
as an additional constraint, the higher-mass stars that do not show offsets from models are 
more likely to hide signals in the low-mass stars that do contain offsets (for example, see 
th e RHS panels of Fig ure [TB]) . It is thus a misleading claim to pronounce the conclusions 



m 



Berger et al.l (120061 ) void, since they are analysin g the wrong mass range (the full span 



of K- and M-dwarfs, as opposed to just M-dwarfs in IBerger et al.ll2006l ). in addition to not 
bringing any r iew rn easu rements for interpreta tion. Lastly, as opposed to the methods in both 
Berger et al.l ( l2006l ) and IDemory et al.l (120091 ) where a solar metallicity model is assumed, we 
compare the observations to models that are specifically tailored to each star's composition. 



5.6. Summary of Global [Luminosity, Temperature, Radius, Mass] Relations 



We find that the empirical relations we derive relating the global properties of: lumi- 
nosity - temperature, luminosity - radius, and temperature - radius show no dependence 
on the metallicity of a star. This conclusion, while unexpected, challenges the Vogt-Russell 
theorem, which states that the properties of a s tar should exhibit a unique solution for a 
star of given mass and chemical composition (see lWeiss et al.ll2004j . and references therein). 
Generally speaking, when regarding relations to global properties of stars, models predict 
that the resulting scatter we observe is caused from varying the stellar metallicity. However, 
the majority of metallicities used in this work are from a uniform source, and should only 
have systematic offsets with respect to each other. So although we do not see this trend for 
the scatter in our data to be caused by metalhcity (e.g., see Figure [13]), we propose that it 
could be a residual effect from underestimated errors in metallicity measurements for these 



^While [Pemorv et all |2009t ) report on the diameters of low-mass stars, only 2 of the stars they observe 
are M-dwarfs, and are repeated measurements of stars with diameters already existing in the literature. 
Furthermore, these points are not averaged, but counted twice in their analysis. 

Pemorv et al. ( 20091 ) quote a 5% error in the derived mass of the star, and thus errors on the theoretical 
radius, and are likely underestimates (see the errors on their figure 9). 
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types of stars (average metallicity error is ~ 0.2 dex). It is important to note however, that 
we are only reviewing this lack of dependence with the composition on the stellar physical 
parameters with respect to iron abundances. Alternate elemental abundances such as Helium 
that are not accounted for in the models can have an effect on a star's physical properties. 
However, our data not only shows no correlations with differences in iron abundances, but 
also suggests that models over-predict the dependence on the stellar global parameters with 
metallicity as illustrated in the pattern of the plotted residuals in the observations versus 
models of Figures [TT] and [T21 

With our data, we are able to develop an empirical method of quantifying the effect 
a star's metallicity has on the observed color index. For instance, our analysis shows that 
metallicity does not effect the stellar radius when relating the global properties of stars 
(luminosity - radius, temperature - radius), while on the other hand, the analysis of the 
color - metallicity - radius relations in Section I^l2] (Figure [8|l . identifies a substantial spread 
in radii for a given color index, noting a dependence on metallicity. This leads us to deduce 
that the metallicity only directly affects the observed color index of a star. The effect can be 
visualized in the color - metallicity - radius relations (Figure |H]) as a shift in the horizontal 
direction for our solutions with different metallicity, as opposed to the vertical direction, 
where at a fixed color index, the radii changes with respect to metallicity. We can quantify 
this shift by defining a zero metallicity color index: we find that for a differential metallicity 
of 0.25 dex, the {B - V), {V - R), {V - I), {V - J), {V - H), and {V - K) color indices 
change by ~ 0.05, 0.15, 0.2, 0.3, 0.3, and 0.3 magnitudes, respectively. This argument, while 
described with respect to our observed radii, also holds using the stellar luminosities as a 
proxy. 



In this work, the masses for sing. 



e stars are derived indirectly via the mass- luminosity 



relation (MLR) in lHenrv fc McCarthvl (Il993r) (Section l3.2p . This relation, as well as others in 
the literature (e.g. IPelfosse et aL 2000), do not incorporate metallicity in their formulation 
and analysis. Therefore, the application of a MLR may produce spurious results for an 
object with a different composition than the objects used to develop the MLRs. To alleviate 
this concern, we used the i^-band form of the MLR, which is thought to be insensitive to 
metallicity: a fortunate by-product to a proper balancing act o f equalizing luminosi ties and 



flux redistributions in the infrared with changing metallicity (see lDelfosse et al.ll2000l for their 
discussion and empirical verification). However, unlike luminosity-temperature, luminosity- 
radius, and temperature-radius planes, we are not able to confirm or deny any metallicity 
dependence on the mass-radius plane with our data because the mass errors are too large 
for the single stars in our sample (Section 15. 5[ Figure [T^ . 



Within Sections 14.11 14.21 and 14.31 we use a multi- variable function to relate the observed 
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colors and metallicities to the global properties of temperature, radius, and luminosity. We 
found that inclusion of metallicity as a variable was essential to properly relate the astro- 
physical quantities (Teff, -R, and L) to the observed colors of our stars. Contrary to the color 
solutions, we find that linking the global properties (temperature, radius, luminosity, mass) 
of a star (15.11 15. 2[ 15. 3[ and 15. 4p do not appear to be sensitive to the stellar metallicity within 
the observational errors used in this work 0. In other words, our data show that the metal- 
licity only affects the observed color-index of a star, and thus using the metallicity-dependent 
transformations in order to convert colors into the stellar properties of temperatures, radii, 
and luminosities is essential. However, the analysis clearly shows that for the range of stars 
that we observe, the metallicity does not impact the physical properties of a star in a way 
that we are able to measure: throwing a bucketful of metals in a star does not make it expand 
in size or cool its surface temperature, it simply morphs the observed color index. 



6. A Need to Validate Binary Star Temperatures 

Typically, when making comparisons of single and binary star properties, the community 
has viewed the data on the mass-radius plane, as we present in Section 15.51 While it is true 
that EEs represent the highest precision in, and largest volume of, the available data for 
stellar masses and radii, it is still not as fundamental when adding single stars to the analysis 
because single stars must rely on mass-luminosity relations to derive their masses. However, 
the data set we introduce here turns the table to the single stars as the most voluminous 
data set of precise measurements to K- and M-dwarf properties. Because of this, we compare 
the fundamental parameters of both the single and binary stars on the temperature-radius 
plane, as they are the most directly determined measurements for single stars. 

While individual component masses and radii can be determined for double-lined EB 
systems, only the ratio of effective temperatures (and thus luminosities), may be calculated. 
Solutions for estimating EB temperatures are generally constrained by photometric calibra - 



tions, where the results often vary due to differences within analyses ([Torres et al.ll2010l ). 
Furthermore, the empirical data to calibrate the temperature scales of late-type stars have 
only recently become available, making any absolute zero point in temperature scales sub- 
ject to uncertainties. The basic technique of EB star temperatures relies on an assumed 
temperature of the primary star that is then translated to the properties of the secondary 
component. The calculation of the EB stellar luminosities are then based on the Stephan- 



^^Since there is a physical coupling of the temperature, radius, and luminosity through the Stephan- 
Boltzmann equation, we expect this connection to their properties to exhibit similar attributes. 
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Boltzmann equation and interstellar extinction (if present). Additionally, a measure of the 
accurate stellar abundances of the EB system is typically unreachable due to the complex 
blending of spectral features, which are dependent on the (unconstrained) individual tem- 



Andersen 


(1991 


) and 


Torres et al. 


(2010) 



Using the same collection of data for binary stars as described in Section 15. 5[ we regard 
the differences on the temperature-radius plane for single stars and binary stars in Figure [T7] 
0. Although there is considerable scatter, the binary star data are clearly separated from 
the single stars on the temperature-radius plane. It appears that either the EB components 
tend to have larger radii (up to ~ 50% for radii < 0.6 M©) compared to a single star of the 
same temperature, or the EB components have cooler surface temperatures (up to ~ 500 K) 
compared to single stars of the same radius. A combination of the two effects could also 
explain the observed disagreement. 

Since our conclusions in the previous section show that mass-radius relationship shows 
minimal differences between single and binary radii (Section 15. 4[ Figure USD, we may infer 
that the systematic differences in Figure [T7] arise from an offset in temperature, not radius, 
between single stars and binary stars. Given an approximate 1:1 relation between mass and 
radius, the observed offset is also present for the data if viewed on the temperature-mass 
plane, where for a given mass, the temperatures of EBs are lower than those of single stars 
by several hundred Kelvin. 

The interpretation of the discrepancy between single and binary star temperatures to be 
a consequence of the diverse physical nature of the two different sets of stars is premature, as 
we must first demonstrate that the observed effect is genuine. As such, even without knowl- 
edge of how the binary stars temperatures are derived in Figure [T71 (technique or calibration), 
it is straightforward to show that depending on what computational method is adopted to 
determine Teff, different solutions for Teff will be obtained. For example. Figure [T71 shows 
the temperature-radius curve derived for single stars (Section |5.3[ Equation [9]) , the 5 Gyr 
solar metallicity Dartmouth model curve, and the tabulated radii and temperatu re values 
for main - sequ e nce st ars in Allen's Astrophysical Quantities (Allen's AQ: ICoxll2000l ). We find 
that the ICoxl (120001 ) temperature values generally agree with the temperatures for the EB 
systems, having temperatures ~ 200 to 300 K cooler than what is predicted by both the 
empirical relation (Equation |9]) and the Dartmouth model. 

Once the method to determine EB temperatures is revisited, we may then proceed 
to uncover if there is any residual discrepancy seen in the temperatures of single and bi- 



^^Note that not all binary stars used in Section [53] have published values for the temperatures. 
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nary stars. Any detection of an offset at this point would likely be of astrophysical origin. 
For inst ance, both highe r mag netic activity and star-spots have be en previously postulated 



by, e.g., [Torres &: RibasI ( 120021 ) and iLopez-Morales fc Ribad (120051 ) . and would result in an 



overall lower effective temperature for a more active star. Based on the normalized X-ra y 
emission Lx/Lbol results from Table [6l and column 5 in table 1 from iLopez- Morales! ( 120071 ). 
the stars in binaries are 10^ — 10^ times more active than the single star counterparts. Such a 
scenario would implicate the elevated activity levels and increase of surface inhomogeneities 
(spots) of EBs playing an important role in the difference of effective temperatures in single 
and binary stars. 



7. Conclusion 

In this work, we present 21 interferometrically determined radii of nearby K and M 
dwarfs. This survey nearly doubles the number of directly determined radii for single stars 
in this range of spectral types. With our measurements of angular diameter and bolometric 
flux, we are able to empirically determine stellar linear radii, luminosities, and effective 
temperatures of the entire sample. Additional properties of the stellar mass and the Lx / Lbol 
activity indicator are computed for each source. We develop empirical relations of the broad- 
band color indices as a functions of temperatures, radii, and luminosities of the objects. We 
find that adding metallicity as an additional parameter is necessary to properly constrain 
our fits, where the data show a notional effect on the color indices due to metallicity among 
the K-type stars, becoming very strong by spectral type ~ MO and later. These relations 
are accurate to ~2% (temperature), ~5% (radius), and ~4% (luminosity). 

Solutions for relating the global stellar parameters (i.e., temperature-radius, temperature- 
luminosity, and mass-radius), on the other hand, are found to be independent of metallicity. 
For example, a solar metallicity star with an observed radius of 0.6 R© will have the same 
mass, temperature, and luminosity as a metal-poor counterpart. This allows us to empiri- 
cally quantify the influence of metallicity on the color index, which appears to impose the 
strongest dependence towards the K- to M- star transition. 

We show that models predict a drop in luminosity at 0.1 Lq and radius at 0.6 R0 at colors 
several tenths of a magnitude too blue than what we observe (Figure[9]). It is unclear whether 
the cause of this is within the model physi cs itself, or whether it lies within the errors of color 



tables used in the color transformations (jPotter et al.l 120081 ). We are able to bypass errors 



introduced from the color tables by looking directly at the empirically measured luminosity- 
temperature and luminosity-radius data. In this, we find that evolutionary models over- 
predict the temperatures for stars with temperatures < 5000 K by ~ 3%, and under-predict 
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the radii for stars with radii < 0.7 by ~ 5%. This follows suit with similar conclusions 
in the literature that models have a tendency to underestimate the radius and overestimate 



the t emperature in order for the luminosity to come out right (jPopperl 119971 : [Morales et al. 



20081 ) ■ Our data also suggest that for the range of stars observed, the influence of metallicity 



on the global parameters appears to be over accounted for in models (Figure [TT] and [T2|) . 

The data set examined in conjunction with a collection of EBs, shows that the radii 
for single stars are consistent within errors to those in eclipsing systems for a given mass, 
both essentially following a 1:1 relation between mass and radius for stars with masses 
< 0.9 Mq and radii < 0.9 R©. On the mass-radius plane, we are able to confirm the 
previously documented offset between the directly measured radii and the ones based on 
theory observed in late- type binary stars at or below the full convection limit of ~ 0.4 Mq. 
We evaluate this discrepancy as a function of mass, metallicity, and activity, revealing a mild 
correlation to the radius offset to the metallicity of a star, where stars with [Fe/H] < —0.35 
all have radii larger than the model predictions (Figure [T6l) . 

Comparing single and binary star properties on the temperature - radius plane, we find 
that for given stellar radius, binary star effective temperature estimates in the literature are 
systematically lower than the interferometrically determined temperatures of single stars by 
~ 200 — 300 K (Figure [T7j) . The major extent of this discrepancy is likely to be caused by 
differences in measuring techniques as well as insufficiently accurate temperature calibrations. 
However, we emphasize that this analysis alone can not exclude the role stellar activity 
may have on stellar temperature and radius. We strongly encourage the validation of EB 
temperatures in order determine to what extent higher activity rates may influence the 
physical properties of binary stars. 

One practical application of measuring these stellar fundamental properties for stars in 
the solar neighborhood is the construction of an empirically determined HR diagram with 
astrophysical, as opposed to observable, axes. This is shown in Figure [181 which comprises 
the 33 stars with diameters presented in this paper as well as those from the literature 
(Table [6]), and the 41 early- type counterparts ( A, F, and G-type dwa rfs) that were examined 



in an analogous way in Paper I of this series (iBoyajian et al.ll2012l ). for a total of 74 stars 



li 



Such a well-populated HR diagram for late-type stars in the theoretical plane, but based 
only on empirical data, has not existed to-date due to the aforementioned paucity in the 
astronomical literature of precise, directly determined radii. 

At the time of writing, the only two stars fro m our compiled sample tha t host exoplanets 



in orbits around them are G J 436 and G J 581 (see von Braun et al. 2012 and von Braun et al. 



20 111 , and references therein). However, ongoing radial velocity and transit exoplanet surveys 



with increasing sensitivity are continuously producing new planet discoveriea^^l. We therefore 
calculate the boundaries of the system habitable zones (HZs) for the stars in Table [6] based 
on our measured stellar astrophysical properties. The HZ boundaries are given in Table [T^ 
where -Rmner and -Router co rrcspoud to th e inner and outer radius of the HZ, respectively, 



following the formalism in iJones fc SleepI ( l2010l ). The HZs for the stars in Paper I are also 
defined in this Table. 

CHARA is currently the only interferometer with the sensitivity and resolution to enable 
this kind of work to this kind of precision. We are thus actively continuing our survey for 
stellar angular diameters to populate the rarefied parameter space of fully characterized 
K- and M- dwarfs in our solar neighborhood. The addition of new data will improve the 
empirical relations we develop in this paper, delivering the foundation needed to facilitate 
fundamental understanding of stellar astrophysics. 

Finally, the ongoing surveys to detect extrasolar planets are making constant progress in 
sensitivity and sophistication, and many of them are concentrating on late-type stars in the 
solar neighborhood because 1) the systems are bright and can thus be studied in great detail, 
and 2) the low stellar mass puts the habitable zone closer to the parent star, making the 
detection of exoplanets in the habitable zones at correspondingly lower periods as compared 
to, e.g., the Earth - Sun system, much more straightforward. Knowledge of the respective 
stellar radius, effective temperature, and luminosity characterizes the radiation environment 
of the orbiting exoplanets and thus gives insight into the first of the crucial conditions that 
must be fullfilled if the planets were to be able to host liquid water on their surfaces. 
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Table 1. Observation Log 





UT 






# of 


Calibrator 


Object 


Date 


Baseline 


Filter 


Brackets 


HD 


GJ 15A 


2008/09 / 16 


El /Wl 


k' 


Q 


6920 




2008 /09 /17 


S 1 /El 


k' 


12 






2008/10/23 


El/Wl 


k' 


3 


905 




2009/12/01 


E1/W2 


K' 


5 


905 




2009 /12/02 


S2/E2 


K' 


5 


3765, 6920 




2010/10/10 


S1/E2 




4 


905 




2011/08/21 


Sl/Wl 




5 


905 


GJ 33 


2009 /08/20 


Sl/El 


k' 


18 


2344, 2454, 6288 




2010 /09 /19 


Sl/El 




5 


4928 


GJ 105 


20t)9 /08 /21 


Sl/El 


k' 


4 


16647, 16970 




2010/09 /16 


El/Wl 


H 


6 


16970 




2010/09/18 


El/Wl 


H 


6 


16970 


nil flfi A 




El/Wl 


k' 


3 


22879 




2008/10/24 


El/Wl 


k' 


6 


22879, 29391 




2009/12/01 


E1/W2 


K' 


5 


22879, 29391 


GJ 205 


2008/10 /23 


El/Wl 


k' 


2 


38858 




2009/12/01 


E1/W2 


k' 


8 


33608, 38858 




2010/09/19 


Sl/El 


H 


12 


33256, 38858 


GJ 338A 


2008/10/22 


Sl/El 


k' 


7 


69548, 80290 




2008/10/24 


El/Wl 


K' 


7 


69548, 80290 




2009/12/01 


E1/W2 


K' 


6 


69548, 80290 




2009/12/02 


S2/E2 


K' 


2 


80290 




2t)09 /12/t)3 


S2/E2 


k' 


2 


69548 


GJ 338B 


2008/10 /22 


Sl/El 


k' 


5 


69548 80290 






El/Wl 


k' 








2009/12/01 


E1/W2 


k' 


6 


69548, 80290 




2009/12/02 


S2/E2 


k' 


1 


80290 




2009 /I 2/03 


S2/E2 


k' 


2 


69548 


GJ 380 


2009/11/22 


Sl/El 


k' 


15 


85795 90508 




2010/11/10 


El/Wl 


H 


3 


89021, 90839 


GJ 411 


2009 /12/02 


S2/E2 


k' 


6 


90277, 99984 




2009 /I 2/03 


S2/E2 


k' 


2 


90277, 99984 




2009 /I 2/04 


S2/E2 


k' 


3 


90840 


GJ 412A 


2009/05/28 


El/Wl 


k' 


5 


99984 




2009/12/02 


S2/E2 


k' 


9 


99984, 103799 




2009/12/03 


S2/E2 


k' 


13 


90277, 99984, 103799 




2009 /I 2/04 


S2/E2 


k' 


11 


90840 99984 103799 


GJ 526 


2006/04/19 


Sl/El 


k' 


11 


119550, 120066 




2009/05/27 


El/Wl 


k' 


2 


121560 




20i()/()(i/:-io 


Sl/El 


H 


4 


119550 




2010/07/01 


Sl/El 


H 


3 


119560 


GJ 631 


2009/05/27 


El/Wl 


k' 


6 


147449 




2009/05/28 


El/Wl 


K' 


7 


150177 




2009/06/18 


Sl/El 


K' 


2 


147449 




2010/04/04 


Sl/El 


H 


5 


147449 


GJ 687 


2004/06/29 


Sl/El 


k' 


1 


161614 




2004/05/27 


Sl/El 


k' 


7 


164633, 151514 




2004/05/31 


Sl/El 


k' 


2 


164633 




2010/06/27 


El/Wl 


H 


6 


168633 




2010/06/28 


El/Wl 


H 


2 


161614 




2010/06/29 


El/Wl 


H 


4 


168633 


GJ 699 


2009/05/28 


El/Wl 


K' 


1 


164363 




2009/08 21 


Sl/El 


K' 


2 


161868 




2010/06/30 


Sl/El 


H 


6 


164363, 171834 
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Table 1 — Continued 





UT 






# of 


Calibrator 


Object 


Date 


Baseline 


Filter 


Brackets 


HD 




2010/07/01 


Sl/El 


H 


4 


164353 


GJ 702A 


2009/05/26 


El/Wl 


k' 


3 


164353 




2009/05/27 


El/Wl 


k' 


2 


164353 




2009/08/21 


Sl/El 


K' 


7 


161868 




2010/07/01 


Sl/Bl 


k' 


8 


164363 


GJ 702B 


2009/05/27 


El/Wl 


k' 


1 


164363 




2009/08/21 


Sl/Bl 


k' 


1 


161868 


GJ 725A 


2008/10/22 


Sl/El 


k' 


3 


168161, 186760 




2009/05/26 


El/Wl 


k' 


3 


168161 




2009/05/27 


El/Wl 


k' 


13 


168161, 173920 


GJ 726B 


2008/10/22 


Sl/El 


K' 


5 


168151, 186760 




2009/05/26 


El/Wl 


k' 


2 


168151 




2009/05/27 


El/Wl 


k' 


14 


168151, 173920 


GJ 809 


2009/11/20 


Sl/El 


K' 


6 


197373 




2009/11/21 


Sl/Bl 


k' 


4 


197373 




2009/11/22 


Sl/Bl 


k' 


4 


197960 




2010/09/16 


Bl/Wl 


H 


10 


197373, 197950 




2010/09/20 


Sl/Bl 


H 


6 


197373, 197950 


GJ 880 


2009/08/21 


Sl/Bl 


k' 


10 


217813, 218261 




2009/12/02 


S2/B2 


k' 


7 


217813, 218261 




2010/07/26 


Sl/Bl 


k' 


10 


217813, 218261 




2010/07/27 


Sl/Bl 


H 


9 


217813, 218235 




2010/07/28 


Sl/Bl 


H 


9 


218236 




2011/08/18 


Sl/Bl 


H 


6 


217813, 218261 


GJ 892 


2008/09/15 


S2/E1 


K' 


10 


219623, 221354 




2008/09/16 


El/Wl 


K' 


9 


219623, 221354 




2008/09/17 


Sl/Bl 


K' 


8 


219623, 221354 




2009/11/22 


Sl/Bl 


K' 


7 


219623, 221354 



-43 - 



Table 2. Calibrators Observed 





Calibrator 


RA 




DEC 


Spectral 


V 


K 


^SED i o- 


Target (s) 


HD 


(hh mm ss) 


(dd mm ss) 


Type 


(mag) 


(mag) 


(mas) 


HD 


905 


00 


13 


31 


41 


02 


07 


FO IV 


5.72 


4.84 


0.397 ± 0.020 


GJ 15A 


2344 


00 


27 


20 


02 


48 


51 


G4 III 


7.19 


4.93 


0.422 ± 0.026 


GJ 33 


2454 


00 


28 


20 


10 


11 


23 


F6 V 


6.04 


4.94, 


0.391 ± 0.014 


GJ 33 


3765 


00 


40 


49 


40 


11 


14 


K2 V 


7.36 


5.16 


0.413 ± 0.035 


GJ 15A 


6288 


01 


03 


49 


01 


22 


01 


Fl V 


6.08 


5.48 


0.290 ± 0.016 


GJ 33 


6920 


01 


10 


19 


42 


04 


53 


F8 V 


5.68 


4.46 


0.508 ± 0.015 


GJ 15A 


16647 


02 


40 


16 


06 


06 


43 


F3 V 


6.26 


5.24 


0.339 ± 0.012 


GJ 105 


16970 


02 


43 


18 


03 


14 


09 


A2 V 


3.47 


3.08 


0.770 ± 0.068 


GJ 105 


22879 


03 


40 


22 


—03 


13 


01 


F9 V 


6.74 


5.18 


0.342 ± 0.021 


GJ 166A 


29391 


04 


37 


36 


—02 


28 


25 


FO V 


5.22 


4.54 


0.432 ± 0.025 


GJ 166A 


33608 


05 


11 


19 


— 02 


29 


27 


F5 V 


5.90 


4.82 


0.415 ± 0.014 


GJ 205 


38858 


05 


48 


35 


—04 


05 


41 


G4 V 


5.97 


4.41 


0.526 ± 0.029 


GJ 205 


69548 


08 


20 


26 


57 


44 


36 


F4 V 


5.89 


4.90 


0.401 ± 0.016 


GJ 338A, GJ 338B 


80290 


09 


20 


44 


51 


15 


58 


F3 V 


6.15 


4.97 


0.385 ± 0.016 


GJ 338A, GJ 338B 


85795 


09 


55 


43 


49 


49 


11 


A3 III 


5.28 


5.08 


0.329 ± 0.031 


GJ 380 




10 


17 


06 


42 


54 


52 








n Tfii -1- n cifiA 

U. IXJI i U.UD41 


GJ 380 


90277 


10 


25 


55 


33 


47 


46 


FO V 


4.73 


3.97 


0.599 ± 0.034 


GJ 411, GJ 412A 


90508 


10 


28 


04 


48 


47 


()(i 


Gl V 


6.44 


4.87 


0.442 ± 0.022 


GJ 380 


90839 


10 


30 


38 


55 


58 


50 


F8 V 


4.83 


3.64 


0.731 ± 0.025 


GJ 380 


90840 


10 


30 


06 


38 


55 


30 


A4 V 


5.79 


5.48 


0.283 ± 0.011 


GJ 411, GJ 412A 


99984 


11 


30 


31 


43 


10 


24 


F4 V 


5.95 


4.59 


0.483 ± 0.020 


GJ 411, GJ 412A 


103799 


11 


57 


15 


40 


20 


37 


F6 V 


6.61 


5.34 


0.343 ± 0.013 


GJ 412A 


119550 


13 


43 


36 


14 


21 


56 


G2 V 


6.94 


5.33 


0.345 ± 0.011 


GJ 526 


120066 


13 


46 


57 


06 


21 


01 


GO V 


6.30 


4.85 


0.428 ± 0.013 


GJ 526 


121560 


13 


55 


50 


14 


03 


23 


F6 V 


6.10 


4.84 


0.422 ± 0.018 


GJ 526 


147449 


16 


22 


04 


01 


01 


44 


FO V 


4.82 


4.09 


0.564 ± 0.022 


GJ 631 


150177 


16 


39 


39 


-09 


33 


16 


FO V 


6.34 


4.98 


0.391 ± 0.019 


GJ 631 


151541 


16 


42 


39 


68 


06 


08 


Kl V 


7.56 


5.69 


0.307 ± 0.023 


GJ 687 


154633 


17 


02 


16 


64 


36 


03 


G5 V 


6.10 


3.89 


0.761 ± 0.057 


GJ 687 


158633 


17 


25 


00 


67 


18 


24 


KO V 


6.43 


4.52 


0.542 ± 0.043 


GJ 687 


161868 


17 


47 


54 


02 


42 


26 


AO V 


3.75 


3.62 


0.638 ± 0.031 


GJ 699, GJ 702A, GJ 702B 


164353 


18 


00 


39 


02 


55 


54 


B5 I 


3.97 


4.00 


0.417 ± 0.030 


GJ 699, GJ 702A, GJ 702B 


168151 


IS 


13 


54 


64 


23 


50 


F5 V 


5.03 


3.94 


0.618 ± 0.023 


GJ 725A, GJ 725B 


170073 


18 


23 


55 


58 


48 


03 


A2 V 


4.99 


4.78 


0.376 ± 0.034 


GJ 725A, GJ 725B 


171834 


18 


36 


39 


06 


40 


IS 


F3 V 


5.44 


4.46 


0.473 ± 0.020 


GJ 699 


173920 


18 


44 


55 


54 


53 


50 


G5 III 


6.26 


4.38 


0.602 ± 0.011 


GJ 725A, GJ 725B 


178207 


19 


04 


55 


53 


23 


48 


AO V 


5.38 


5.41 


0.276 ± 0.023 


GJ 725A, GJ 725B 


186760 


19 


43 


14 


58 


00 


60 


GO V 


6.30 


4.91 


0.431 ± 0.018 


GJ 725A, GJ 725B 


197373 


20 


40 


18 


60 


30 


19 


F6 V 


6.03 


4.94 


0.404 ± 0.013 


GJ 809 


197950 


20 


43 


11 


66 


39 


27 


AS V 


5.60 


5.06 


0.349 ± 0.012 


GJ 809 


217813 


23 


03 


05 


20 


55 


07 


Gl V 


6.66 


5.15 


0.377 ± 0.013 


GJ 880 


218235 


23 


06 


18 


18 


31 


04 


F6 V 


6.13 


5.07 


0.379 ± 0.010 


GJ 880 


218261 


23 


06 


32 


19 


54 


39 


F7 V 


6.30 


5.14 


0.380 ± 0.012 


GJ 880 


219623 


23 


16 


42 


53 


12 


49 


F7 V 


5.60 


4.31 


0.518 ± 0.027 


GJ 892 


221354 


23 


31 


22 


59 


09 


56 


K2 V 


6.74 


4.80 


0.451 ± 0.020 


GJ 892 



-44- 



Table 3. Angular Diameters of K and M Dwarfs 



Star 9^ of Reduced ^UD i o" ^LD i ^LD 

Name Obs. (mas) fi^ (mas) % err 



GJ 


15A 


37 





.93 


0, 


,970 


± 


0. 


,005 


0, 


,398 


1, 


,005 


± 


0. 


,005 


0, 


,5 


GJ 


33 


23 





.86 


0, 


,836 


± 


0. 


,004 


0, 


,425 


0, 


,868 


± 


0. 


,004 


0, 


,4 


GJ 


105 


16 





.76 


0, 


,989 


± 


0. 


,007 


0, 


,443 


1, 


,030 


± 


0. 


,007 


0, 


,7 


GJ 


166A 


14 


1 


.01 


1 


444 


± 


0. 


,006 


0, 


,410 


1, 


,504 


± 


0. 


,006 


0, 


,4 


GJ 


205 


22 


1 


.40 


0, 


,904 


± 


0. 


,003 


0, 


,453 


0, 


,943 


± 


0. 


,004 


0, 


,4 


GJ 


338A 


23 





.74 


0, 


,834 


± 


0. 


014 


0, 


,491 


0, 


,871 


± 


0. 


,015 


1 


,7 


GJ 


338B 


21 





.76 


0, 


,823 


± 


0. 


,016 


0, 


,453 


0, 


,856 


± 


0. 


,016 


1 


,9 


GJ 


380 


18 


1 


.00 


1, 


169 


± 


0. 


,008 


0, 


,491 


1, 


225 


± 


0. 


,009 


0, 


,7 


GJ 


411 


11 


1 


.75 


1, 


,380 


± 


0. 


,013 


0, 


,391 


1, 


432 


± 


0. 


,013 


0, 


,9 


GJ 


412A 


35 


1 


.25 


0, 


739 


± 


0. 


,016 


0, 


,398 


0, 


764 


± 


0. 


017 


2 


,2 


GJ 


526 


29 


2, 


.60 


0, 


,807 


± 


0. 


,013 


0, 


,398 


0, 


,835 


± 


0. 


014 


1 


,6 


GJ 


631 


20 


0, 


.54 


0, 


701 


± 


0. 


Oil 


0, 


,396 


0, 


724 


± 


0. 


Oil 


1 


,6 


GJ 


687 


22 


1 


.63 


0, 


,830 


± 


0. 


,013 


0, 


,391 


0, 


,859 


± 


0. 


014 


1 


,6 


GJ 


699 


13 





.66 


0, 


917 


± 


0. 


,005 


0, 


,408 


0, 


,952 


± 


0. 


,005 


0, 


,6 


GJ 


702A 


20 


1 


.51 


1, 


,460 


± 


0. 


,004 


0, 


,396 


1, 


,515 


± 


0. 


,005 


0, 


,3 


GJ 


702B 


2 





.35 


1, 


169 


± 


0. 


,015 


0, 


,476 


1, 


221 


± 


0. 


,015 


1 


,3 


GJ 


725A 


19 





.58 


0, 


,907 


± 


0. 


,008 


0, 


,391 


0, 


,937 


± 


0. 


,008 


0, 


,9 


GJ 


725B 


21 


1 


.42 


0, 


822 


± 


0. 


,015 


0, 


,408 


0, 


,851 


± 


0. 


,015 


1 


,8 


GJ 


809 


31 


1 


.33 


0, 


,698 


± 


0. 


,008 


0, 


,398 


0, 


722 


± 


0. 


,008 


1 


,1 


GJ 


880 


51 





.94 


0, 


716 


± 


0. 


,004 


0, 


,453 


0, 


744 


± 


0. 


,004 


0, 


,5 


GJ 


892 


34 





.69 


1, 


,063 


± 


0. 


,007 


0, 


,443 


1 


106 


± 


0. 


,007 


0, 


,7 



Table 4. Comparison of Angular Diameters 



star 
Name 




9ld ± <^ 
(mas) 




Reference 


Instrument 


Spectral 
Type 


% err 






GJ 15A 


1.005 


± 0.005 




This work 


CHARA 


Ml. 5 V 


0, 


5 


0, 


,0 




0, 


,988 


± 0.016 




Borser et al. f2006) 


CHARA 




1, 


6 


1 


,0 




1.027 ± 0.059 


van 


Belle & von Braun (2009) 


PTI 




5, 


7"= 


-0, 


,4 


GJ 33 


0.868 ± 0.004 




This work 


CHARA 


K2 V 


0, 


5 


0, 


.0 




0.933 ± 0.064 


van 


Belle & von Braun f 20091 


PTI 




6. 


,9*= 


-1 


.0 


GJ 105 


1.03 ± 0.007 




This work 


CHARA 


K3 V 


0, 


7 


0, 


.0 




0, 


.936 


± 


0, 


.070 




Lane et al. (20011 


PTI 




7, 


5*= 


1 


.3 


GJ 205 


0, 


,943 


± 


0, 


.004 




This work 


CHARA 


Ml. 5 V 


0, 


4 


0, 


.0 




1 


,149 


± 


0, 


,110 




Seeransan et al. (2003) 


VLTI 




9, 


,6'= 


-1 


.9 


GJ 166A 


1 


.504 


± 


0, 


.008 




This work 


CHARA 


Kl Vc 


0, 


4 


0, 


.0 




1 


.437 


± 


0, 


.039 




Demorv et al. (20091 


VLTI 




2, 


7 


1 


.7 


GJ 380 


1 


.225 


± 


0, 


.008 




This work 


CHARA 


K7.0 V 


0, 


7 


0, 


.0 




1 


.155 


± 


0, 


.040*= 




Lane et al. (2001) 


PTI 




3, 


5 


1 


.7 




1 


.238 


± 


0, 


.053 


van 


Belle & von Braun (20091 


PTI 




4, 


,3 


-0, 


.2 


GJ 411 


1 


.432 


± 


0, 


.013 




This work 


CHARA 


M2.0 V 


0, 


,9 


0, 


.0 




1 


.436 


± 


0, 


.030 




Lane et al. (2001) 


PTI 




2, 


1 


-0, 


1 




1 


,439 


± 


0, 


,048 


van 


Belle & von Braun (20091 


PTI 




3, 


,3 


-0, 


.1 


GJ 526 


0, 


,835 


± 


0, 


,014 




This work 


CHARA 


Ml. 5 V 


1, 


7 


0, 


.0 




0, 


.845 


± 


0, 


.057 




Bereer et al. (20061 


CHARA 




6. 


.7" 


-0, 


.2 


GJ 631 


0, 


.724 


± 


0, 


.011 




This work 


CHARA 


KO V 


1. 


,5 


0, 


.0 




0, 


.888 


± 


0, 


.066 


van 


Belle & von Braun (20091 


PTI 




7, 


,4'= 


-2 


.5 


GJ 687 


0, 


.859 


± 


0, 


.014 




This work 


CHARA 


M3.0 V 


1, 


,6 


0, 


.0 




1 


.009 


± 


0, 


.077 




Berser et al. (20061 


CHARA 




7, 


,6'= 


-1 


.9 


GJ 699 


0, 


,952 


± 


0, 


.005 




This work 


CHARA 


M4.0 V 


0, 


5 


0, 


.0 




1 


,004 


± 


0, 


.040 




Lane et al. (20011 


PTI 




4, 


,0 


-1 


.3 


GJ 880 


0, 


,744 


± 


0, 


.004 




This work 


CHARA 


Ml. 5 V 


0, 


5 


0, 


.0 




0, 


,934 


± 


0, 


,059 




Berser et al. (20061 


CHARA 




6, 


,3'= 


-3, 


.2 



^Wc define the combined error as ctq — ("^xhiswork "^Reference ' 

'^Wc caution the reader that the value of cited in lLane et all i200l| ) is less than 9\jjj. 

■^These measurements have errors > 5%. See Scction[^for details. 



Table 5. Object Photometry used in Relations 





B 


V 


R 


I 


J 


H 


K 


star 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


GJ 15A . 


9.63 


8.07 


6.69 


5.53 


4.86 


4.25 


4.02 


GJ 33 . . 


6.64 


5.76 


4.99 


4.52 


4.24 


3.72 


3.61 


GJ 53A . 


5.87 


5.18 


4.55 


4.14 


3.86 


3.39 


3.36 


GJ 75 . . 


6.44 


5.63 


4.99 


4.60 


4.31 


3.88 


3.84 


GJ 105 . 


6.79 


5.82 


4.99 


4.46 


4.07 


3.52 


3.45 


GJ 144 . 


4.62 


3.74 


3.01 


2.55 


2.20 


1.75 


1.65 


GJ 166A 


5.25 


4.43 


3.74 


3.29 


2.95 


2.48 


2.41 


GJ 205 . 


9.44 


7.97 


6.53 


5.39 


4.77 


4.06 


3.86 


GJ 338A 


9.05 


7.64 






4.89 


4.25 


4.09 


GJ 338B 


9.04 


7.70 






4.78 


4.30 


4.15 


GJ 380 . 


7.94 


6.59 


5.36 


4.56 


3.98 


3.32 


3.19 


GJ 411 . 


9.00 


7.49 


5.99 


4.80 


4.13 


3.56 


3.35 


GJ 412A 


10.32 


8.77 


7.29 


6.22 


5.56 


4.95 


4.76 


GJ 436 . 


12.17 


10.65 






6.90 


6.32 


6.07 


GJ 526 . 


9.93 


8.50 


7.06 


5.92 


5.26 


4.64 


4.46 


GJ 551 . 


13.02 


11.05 


8.68 


6.42 


5.33 


4.73 


4.37 


GJ 570A 


6.88 


5.78 


4.85 


4.28 


3.82 


3.27 


3.15 


GJ 581 . 


12.19 


10.58 


8.89 


7.46 


6.68 


6.09 


5.83 


GJ 631 . 


6.55 


5.74 


5.13 


4.74 


4.32 


3.86 


3.83 


GJ 687 . 


10.65 


9.15 






5.37 


4.75 


4.54 


GJ 699 . 


11.27 


9.54 






5.30 


4.77 


4.50 


GJ 702A 


5.06 


4.20 










2.30 


GJ 702B 


7.15 


6.05 










3.10 


GJ 725A 


10.44 


8.90 






5.19 


4.68 


4.46 


GJ 725B 


11.28 


9.69 






5.72 


5.24 


4.98 


GJ 764 . 


5.49 


4.69 


4.04 


3.63 


3.32 


3.04 


2.78 


GJ 809 . 


10.04 


8.58 


7.20 


6.14 


5.52 


4.81 


4.64 


GJ 820A 


6.39 


5.23 


4.20 


3.56 


3.16 


2.61 


2.40 


GJ 820B 


7.38 


6.02 


4.87 


4.07 


3.58 


2.93 


2.73 


GJ 845 . 


5.75 


4.69 


3.81 


3.25 


2.83 


2.30 


2.18 


GJ 880 . 


10.19 


8.68 






5.41 


4.78 


4.58 


GJ 887 . 


8.83 


7.35 






4.20 


3.60 


3.36 


GJ 892 . 


6.57 


5.57 


4.74 


4.21 


3.86 


3.40 


3.23 



R,eference^\//^/ , Reference j/fx 



JohnBon fl9651. Glass (1975 ^ 
JohnsonetaL (1966) , Johnson ct al. (1968) 
Johnsonet^. ( 1966) .Johnson ct al. (1968) 

Johnson et al. (1966 ) .Glass (1975 1 
Johnsoi^_et_aL Q966 ) ,Glass ^975) 
^^^^^o^^^^^. (1966) .Johnson ct al. ^^68) 
Johnson Q965). Mould & Hyland ( 1976 1 
Cowley et al. ( 1967) , Vccdcr (1974) 
Cowley et al. (1967).VcodCT (1974) 
Johnson ^1965), Glass (1975 1 
Johnson [1964), Glass (1975 ) 
Johnson (^965), Glass (1975 ) 
^cis (1993).Cutri ct al. '(2003b) 
"johnson (1965), Glass (1975)' 
Frogel et al. (1972 ), Mould & Hyland (1976) 
Johnson (1965 1, Mould & Hyland (1976)'" 
Johnson (1965 i,M ould & Hyland (1976 ) 
JohnsonetaL (1966 ) jJohnsonetaL (1968) 
^o^so^^^Mo^^ ^53 ) ^e^so^^^^ ^77) 
^ormilliod (1997), Mould & Hyland (1976) 
Eggenbcrgcr ct al. (2008b) .Christen & Drummond [2006) 
Eggen^crgcr ct al. (2008b) .Chris ten & Drummond [2006) 
Johnson & Morgan (195 3) .Vccdcr (1974 ) 
Johnson &: Morgan (1953) j^cdcr [1974 ) 
^^^^^^^^^al. (1966) ,,To^^^^^^^^ ^^66) 

Erro (1971),Pcrsso n et al. (1977 ) 
^^^^^o^^^^^^ ^^^6)^^inson ct al. (1968J 
^^^^^o^^^^^^ ^^66)^^inson ct al. (1968) 
Johnson et al. fl966). M ould fc Hyland [1976) 
Err o (1971), Erro (1971 ) 
Bessel (1990 ) , Mould & Hyland (1976 ) 
JohnsonetaL [1966 ) jJohnsonetaL [1968) 



Note. — 



Refer to Section 13.21 for details. 



Table 6. Fundamental Parameters 



star 


Spectral 


Metallicity 


Metallicity 




R± 






Radius 






L ±a 


^EFF it 




Mass 






Name 


Type 


[Fe/H] 


Reference^ 




(«©) 




Reference+' 


(10-* erg/cm/s^) 






(K) 






ix/^'BOL 


GJ 15At 


Ml. 6 V 


— 0. 


.36 


1 


0, 


.3874 


± 


0. 


.0023 


this work 


5.420 ± 0.044 


0. 


,02173 


± 0.00021 


3663 


± 


11 





,423 


2 


.32E-06 


GJ 33 


K2 V 


— 0. 


.22 


2 


0, 


.6964 


± 


0. 


.0041 


this work 


16.060 ± 0.100 


0. 


,26073 


± 0.00263 


4960 


± 


14 





.763 


8, 


.69E-07 


GJ 105 


K3 V 


-0 


.08 


2 


0, 


.7949 


± 


0. 


,0062 


this work 


16.680 ± 0.103 


0, 


,26790 


it 0.00239 


4662 


± 


17 


0, 


.767 


1 


.70E-06 


GJ 166At 


Kl Vc 


-0 


.24 


2 


0, 


.8061 


± 





0036 


this work 


52.690 ± 0.394 





40782 


it 0.00319 


5143 


it 


14 


0, 


.816 


1 


.03E-05 


GJ 205 


Ml. 5 V 


0. 


.35 


3 


0, 


.5735 


± 


0. 


.0044 


this work 


6.182 ± 0.032 


0, 


,06163 


± 0.00088 


3801 


± 


9 


0, 


.615 


1 


.62E-05 


GJ 338A 


MO.O V 


— 0. 


.18 


3 


0, 


.6773 


± 


0. 


.0131 


this work 


5.885 ± 0.061 


0. 


,06974 


± 0.00213 


3907 


± 


36 





.622 


2 


.07E-06 


GJ 338B 


K7.0 V 


— 0. 


.15 


3 


0, 


.6673 


± 


0. 


.0137 


this work 


5.455 ± 0.037 


0. 


,06465 


± 0.00194 


3867 


± 


37 





.600 


2 


.23E-06 


GJ 380* 


K7.0 V 


-0. 


.16 


2 


0, 


.6416 


± 


0. 


.0048 


this work 


13.860 ± 0.093 


0. 


,10253 


± 0.00088 


4081 


± 


16 


0, 


.660 


6, 


.92E-06 


GJ 41lt 


M2.0 V 


-0, 


.41 


3 


0, 


.3921 


± 


0, 


.0037 


this work 


9.842 ± 0.060 


0. 


,01989 


± 0.00014 


3465 


± 


17 


0, 


.403 


8 


.03E-06 


GJ 412A 


Ml.O V 


-0 


.40 


3 


0, 


.3982 


± 


0. 


0091 


this work 


2.908 ± 0.022 


0, 


,02129 


it 0.00026 


3497 


it 


39 


0, 


.403 


3 


.llE-05 


GJ 526 


Ml. 5 V 


— 


.30 


3 


0, 


.4840 




0. 


0084 


this work 


3.979 ± 0.030 


0, 


,03603 


it 0.00051 


3618 


it 


31 





.520 


5 


.33E-06 


GJ 631 


KO V 





.04 


2 


0, 


.7591 


± 


0. 


0122 


this work 


14.160 ± 0.090 


0, 


41945 


it 0.00422 


5337 


it 


41 





.821 


8 


.92E-06 


GJ 687 


M3.0 V 


— 


.09 


3 


0, 


.4183 


it 


0. 


0070 


this work 


3.332 ± 0.022 





02128 


it 0.00023 


3413 


it 


28 


0, 


.413 




. 30E-06 


GJ 699+ 


M4.0 V 


— 0. 


.39 


3 


0, 


. 1867 


i 


0. 


.0012 


this work 


3.262 ± 0.022 


0, 


,00338 


it 0.00003 


Snt 


it 


10 





a^K 


3 


.yyE-uo 


GJ 702A 


KO Ve 






2 


0, 




_i_ 
It 




.0044 


this work 






0.53 ± 0.02= 




_j_ 


62 








.67E-06 


GJ 702B 


K5 Ve 


n 


m 


2 


0, 


6697 


_j_ 


n 


.0089 


this work 






0.15 ± 0.02"= 




_j_ 


149 


n 


698 




.yuEj-UD 


GJ 725A 


MS.C V 




49 


3 


0, 


3561 


_j_ 




0039 


this work 


3.937 ± 0.004 


0, 


,01531 


it 0.00018 


3407 


_j_ 


15 




318 




. 05E-06 


GJ 725B 


MS. 6 V 


Q 


36 


3 


0, 


3232 


_j_ 


g 


0061 


this work 


2.238 ± 0.013 


0, 


,00871 


it 0.00012 


3104 


_j_ 


28 





235 




. OOI^-UQ 


GJ 809 


MO. 5 


— 


.21 


3 


0, 


.5472 


-1- 


0. 


0067 


this work 


3.224 ± 0.027 


0, 


,04990 


it 0.00062 


3692 


-1- 


22 





.573 




.41E-05 


GJ 880 


Ml. 5 V 





.06 


1 


0, 


.5477 


± 


0. 


0048 


this work 


3.502 ± 0.017 


0, 


05112 


it 0.00074 


3713 


± 


11 





.569 


6 


.66E-06 


GJ 892 


K3 V 


0. 


.07 


2 


0, 


.7784 




0. 


.0053 


this work 


19.850 ± 0.086 


0, 


,26499 


it 0.00152 


4699 


± 


16 





.763 


3 


.73E-07 


— n 

GJ 15AT 


M1.5V 


— 


.36 


1 


0, 


.3790 


-1- 


0. 


.0060 


1 


6.420 ± 0.044 


0. 


,02173 


± 0.00017 


3694 


-1- 


30 





.423 


2 


.32E-06 


GJ 53A 


G6Vp 


— 0. 


.68 


4 


0. 


.7910 


± 


0. 


.0080 


3 


25.790 ± 0.147 


0. 


,45845 


± 0.00260 


6348 


± 


26 





.808 


3 


.73E-07 


GJ 75 


KOV 


0. 


.03 


2 


0. 


.8190 


± 


0. 


.0240 


3 


16.460 ± 0.126 


0. 


,61971 


± 0.00396 


6398 


± 


76 


0, 


.837 


1 


.12E-06 


GJ 144 


K2V 


-0. 


.06 


2 


0. 


.7350 


± 


0. 


0050 


5 


102.100 ± 0.457 


0. 


,32897 


± 0.00147 


6077 


± 


36 


0, 


.781 


1 


.64E-06 


GJ 166A ' 


KIV 


-0. 


.24 


2 





.7700 


± 


0. 


0210 


6 


52.690 it 0.394 


0. 


,40782 


it 0.00304 


5261 


± 


72 


0, 


.818 


1 


.03E-05 


GJ 380+ 


K7V 


-0, 


.16 


2 





.6050 


± 


0. 


0200 


4 


13.860 ± 0.093 


0. 


10253 


± 0.00069 


4203 


it 


73 


0, 


.669 


6 


.92E-06 


GJ 380+ 


K7V 


— 0, 


.16 


2 


0, 


.6490 


± 


0. 


0280 


2 


13.860 it 0.093 


0, 


10253 


± 0.00069 


4060 


± 


87 


0, 


.669 


6 


.92E-06 


GJ 411 + 


M2V 


-0, 


.41 


3 


0, 


.3930 


± 


0. 


0080 


4 


9.842 it 0.060 


0. 


01989 


± 0.00012 


3460 


± 


37 





.405 


8 


.04E-06 


GJ 411 + 


M2V 


-0, 


.41 


3 


0. 


.3950 


± 


0. 


0130 


2 


9.842 ± 0.060 


0. 


,01989 


± 0.00012 


3457 


± 


58 


0, 


.405 


8, 


.04E-06 


GJ 436 


M3V 


0, 


.04 


3 


0. 


.4646 


± 


0. 


.0182 


10 


0.788 ± 0.004 


0. 


,02626 


± 0.00012 


3416 


± 


63 


0, 


.472 


6, 


.84E-06 


GJ 551 


M6.6V 


0, 


.19 


6 


0, 


.1410 


± 


0. 


.0070 


6 


2.961 ± 0.037 


0. 


,00166 


± 0.00002 


3064 


± 


79 


0, 


.118 


2, 


.83E-04 


GJ 570A 


K4V 


0, 


.02 


2 


0, 


.7390 


± 


0. 


0190 


6 


19.040 ± 0.164 


0. 


20232 


± 0.00163 


4607 


± 


68 


0, 


.743 


2, 


.72E-06 


GJ 581 


M2.6V 


-0, 


.10 


3 


0, 


.2990 


± 


0. 


,0100 


9 


0.930 ± 0.006 


0. 


,01130 


± 0.00008 


3442 


± 


64 


0, 


.297 


8, 


.28E-06 


GJ 699+ 


M4Ve 


-0, 


.39 


3 


0, 


.1960 


± 


0, 


,0080 


4 


3.262 it 0.022 


0, 


,00338 


it 0.00002 


3140 


± 


63 


0, 


.150 


3, 


.99E-06 


GJ 764 


KOV 


-0, 


.19 


2 


0, 


.7780 


± 


0, 


,0080 


3 


39.670 ± 0.194 


0, 


,40926 


± 0.00199 


5246 


± 


26 


0, 


.806 


2 


.61E-06 


GJ 820A+ + 


K5V 


-0, 


.19 


2 





.6650 


± 


0, 


,0050 


8 


37.750 ± 0.188 


0, 


14295 


± 0.00071 


4355 


it 


17 


0, 


.680 


5, 


.OOE-06 


GJ 820Att 


K5V 


-0, 


.19 


2 





6100 


± 


0, 


,0180 


2 


37.750 ± 0.188 


0, 


14295 


± 0.00071 


4548 


it 


64 


0, 


.680 


5, 


.OOE-06 


GJ 820B+ + 


K7V 


-0, 


.29 


2 


0, 


.5950 


± 


0. 


,0080 


8 


20.340 ± 0.107 


0. 


,07753 


± 0.00041 


3954 


± 


28 


0, 


.629 


9, 


.28E-06 


GJ 820B+ + 


K7V 


-0, 


.29 


2 


0, 


.6280 


± 


0. 


,0170 


2 


20.340 ± 0.107 


0. 


,07763 


± 0.00041 


3862 


± 


63 


0, 


.629 


9, 


.28E-06 


GJ 845 


K5V 


-0, 


.07 


2 


0, 


.7320 


± 


0, 


,0060 


6 


50.730 ± 0.505 


0, 


,20737 


± 0.00206 


4555 


± 


24 


0, 


.731 


1 


.94E-06 


GJ 887+ + 


Mn.r,\- 


n 


1') 




n 


.i')in 







01 10 




1().<)2() ± n.127 







± ().()()(),12 


:!fil2 






n 






.2,8E,-()(i 


GJ SS7 ' 




— u 


i!) 


2 


u 


4,';!)ij 


± 


0. 


()ii(J 


I) 


i().<)2() ± lJ.i27 


0, 


OliO.ji 


± (J.(J()()42 


:!727 


± 


47 


u 


.'522 




,28E-()(i 



su 



,t,tt 



< R> ±a (Rq) 



GJ ISA 
GJ 166A 
GJ 380 
GJ 411 



0.3863it0.0021 
0.8051±0.0035 
0.6398±0.0046 
0.3924±0.0033 



< Teff > ±cr (K) 



Average from 



3567±11 
5147±14 
4085±14 
3464±16 



this work, 1 
this work, 2 
this work, 2, 4 
this work, 2, 4 



Table 6 — Continued 



star 
Name 


Spectral 
Type 


Metallicity 
[Fe/H] 


Metallicity 
Reference'^ 


R± <T 


Radius ^^BOL i 
Reference^ (10~^ erg/cm/s^) 


L ±a 


TefF i 
(K) 


Mass 


-^x/-^BOL 


GJ 699 
GJ 820A 
GJ 820B 
GJ 887 








0.1869±0.0012 
0.6611±0.0048 
0.6010±0.0072 
0.4712±0.0086 






3222±10 
4361±17 
3932±25 
3676±35 




this work, 4 
2, 8 
2, 8 
6, 7 



"M etallicity references are l UNeves et al.H201lH . 21 lAnderson fc Franci J i201ll 1 . 31 lRoias-Avala et ani2012l 1. 41 lBovaiian et al.H2008l 1. SI lEdvardsson et anil993l 1. and 6) lBonflls et all 
i2005l 1. See Section [X2l for details. 

^In terfe rometric references f or measured radii of stars 1) iBerger et al.l ^2006^. 2 ) Ivan Belle fc von BraunI i2009l) . 3) | B ovaiian et al.l 120061) . 4) iLane et all i200ll) . 5) Idi Folco et al.l 
120071 ) , 6) iDemorv et al.l i2009l) , 7) ISegransan et al.l i2003D , SI lKervella et al.l i200el 'l , 9 Hvon Braun et al.l i201ll 'l. 10 1 Ivon Braun et al.l i2012l ). 

'^The photometry for these two sources is too blended to perform this measurement for this work. This value if from fEggenberger et al ] i200Sal1 . 



Note. — The top portion are new measurements made in this work. The middle portion lists measurements of stellar radii found in the literature, with precision of better than 5%. 
Stars with multiple measurements are marked with a ^ or a Stars are marked with a ^ if the mean is includes a measurement from this work, and a if the mean is from two 

literature measurements. The bottom portion of the table lists the stars with multiple measurements, and the weighted mean for their radii and temperatures, all other parameters 
remain unaffected when combining the multiple sources for measured radii. All bolometric flux, luminosity, temperature, mass, and i^x/^BOL values are computed/measured in this 
work. Refer to Sections [3l fsTsl and[^for details. 



Table 7. G-type stars from 



Boyajian et al 



(120121 ) 



star 


B 


V 


R 


/ 


J 


H 


K 


^EFF ± o- 


Metallicity 




HD 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(K) 


[Fe/H] 


R-eference^^/jj , Reference j/j^^ 


19373 


4.65 


4.05 


3.52 


3.23 


3.06 


2.73 


2.69 


5915±29 


0.09 


Johnson (1964); Johnson^^^^ (_1968) 


34411 


5.33 


4.71 


4.18 


3.86 


3.62 


3.33 


3.28 


5 749 ±48 


0.05 


^Johnson (1964); Johnson^^^ (_1968) 


39587 


5.00 


4.41 


3.90 


3.59 


3.34 


3.04 


2.97 


5961±36 


-0.16 


^Johnson (1964); Johnson^^^ (_1968) 


82885 


6.18 


5.41 


4.79 


4.42 


4.14 


3.77 


3.70 


5434±45 


0.06 


Johnson (1964); Ducati (2002) 


101501 


6.08 


5.34 


4.73 


4.37 


4.02 


3.61 


3.60 


5270±32 


-0.12 


Johnson (1964); Johnson et al. (1968) 


109358 


4.86 


4.27 


3.73 


3.42 


3.20 


2.88 


2.83 


5653±72 


-0.30 


Johnson (1964); Ducati (2002) 


114710 


4.84 


4.26 


3.77 


3.47 


3.22 


2.95 


2.89 


5936±33 


-0.06 


Johnson (1964); Johnson et al. (1968) 



Note. — Photometry, temperatures, and metallicities of G dwarfs in lBoyaiian et al ] )2012D . Refer to S ection[4. l|for details. 



Table 8. Solutions to Teff Relations 



Coefficient 


(S - V) - [Fe/H] 


{V - R) - [Fe/H] 


(V - /) - [Fe/H] 


(V - J) - [Fe/H] 


(V -H) - [Fe/H] 


iy - K) - [Fe/H] 


ao 


8010 ± 52 


7646 ± 42 


7325 ± 33 


7308 ± 26 


7641 ± 33 


7643 ± 32 


ai 


-4095 ± 94 


-4295 ± 80 


-2262 ± 35 


-1775 ± 21 


-1611 ± 22 


-1523 ± 21 


12 


819 ± 41 


1058 ± 37 


313 ± 9 


198 ± 4 


151 ± 4 


134 ± 3 


13 


133 ± 57 


304 ± 66 


-15 ± 34 


71 ± 20 


177 ± 19 


137 ± 17 


a4 


39 ± 85 


-77 ± 92 


393 ± 84 


100 ± 64 


-319 ± 73 


-202 ± 72 


15 


-362 ± 57 


102 ± 78 


733 ± 77 


317 ± 54 


185 ± 55 


157 ± 56 


n 

Range 

Median dT (K) 


GKM=40; KM=33 

0.7 - 1.7 
GKM=70; KM=68 


GKM=29; KM=22 

0.6 - 1.7 
GKM=63; KM=62 


GKM=29; KM=22 

1.0 - 3.0 
GKM=73; KM=68 


GKM=38; KM=31 

1.2 - 4.0 
GKM=75; KM=70 


GKM=38; KM=31 

1.5 - 4.5 
GKM=49; KM=48 


GKM=40; KM=33 

1.5 - 5.0 
GKM=49; KM=43 



Note. — The values , a\, a2, 0,3, 0,4, ^5 are the coefficients to the color-metallicity-tcmperaturc relations in the form of Equation [S] The bottom three 
rows show n, the number of data points used for the combined G, K, and M star fit as well as just the K and M star fit, Range, the range in color index 
where each relation holds true, and Median dT, the median absolute deviation for the combined G, K, and M star fit as well as just the K and M star fit. 
See Figurc[6]for a graphical representation of the solutions, and refer to Section |4. 1| for discussion. 



Table 9. Solutions to Radius Relations 



Coefficient 


(S - V) - [Fe/H] 


(V - R) - [Fe/H] 


(V - 7) - [Fe/H] 


(V - J) - [Fe/H] 


(V - H) - [Fe/H] 


(V - K) - [Fe/H] 


10 


0.3830 ± 0.0174 


0.9122 ± 0.0218 


0.8632 ± 0.0180 


0.9647 ± 0.0082 


0.8843 ± 0.0107 


0.9191 ± 0.0097 


ai 


0.9907 ± 0.0297 


-0.0936 ± 0.0420 


-0.0198 ± 0.0199 


-0.1018 ± 0.0058 


-0.0106 ± 0.0067 


-0.0230 ± 0.0060 


12 


-0.6038 ± 0.0124 


-0.1274 ± 0.0196 


-0.0472 ± 0.0052 


-0.0125 ± 0.0011 


-0.0222 ± 0.0011 


-0.0191 ± 0.0009 


^3 


0.1470 ± 0.0236 


0.3682 ± 0.0250 


0.1566 ± 0.0129 


0.0786 ± 0.0073 


0.1064 ± 0.0068 


0.0789 ± 0.0058 


a4 


-0.0451 ± 0.0340 


-0.2662 ± 0.0356 


-0.1531 ± 0.0328 


-0.0179 ± 0.0240 


-0.1533 ± 0.0271 


-0.0643 ± 0.0239 


a5 


-0.1340 ± 0.0222 


-0.1652 ± 0.0279 


-0.0321 ± 0.0277 


0.0286 ± 0.0187 


-0.0127 ± 0.0193 


0.0089 ± 0.0189 


n 


31 


21 


21 


30 


30 


32 


Range 


0.8 - 1.6 


0.6 - 1.6 


1.0 - 3.2 


1.3 - 4.3 


1.6 - 4.8 


1.7 - 5.0 


Median dR (Rq) 


0.030 


0.041 


0.026 


0.036 


0.027 


0.027 



O 



Note. — The values , , a2 . a3 , a4 , 05 are the coefficients to the color-metallicity-radius relations in the form of Equation[^ The bottom three rows show 
n, the number of data points used in the fit, Range, the range in color index where each relation holds true, and Median di?, the median absolute deviation of 
the fit. See Figure[8]for a graphical representation of the solutions, and refer to Section 14.21 for details. 



Table 10. Solutions to Luminosity Relations 



Coefficient (S - V) - [FeH] (V ~ B.) - [FeH] (V - I) - [FeH] (V - J) - [FeH] (_V - H) - [FeH] (V ~ K) - [FeH] 



ao -1.2895 ± 0.0379 0.0132 ± 0.0257 0.1838 ± 0.0208 0.1571 ± 0.0157 -0.0851 ± 0.0208 0.0613 ± 0.0196 

ai 2.5206 ± 0.0667 -0.3282 ± 0.0498 -0.4935 ± 0.0232 -0.3525 ± 0.0129 -0.0494 ± 0.0143 -0.1329 ± 0.0129 

as -1.7633 ± 0.0282 -0.4462 ± 0.0225 -0.0476 ± 0.0059 -0.0364 ± 0.0025 -0.0726 ± 0.0023 -0.0535 ± 0.0020 

a-i 0.4953 ± 0.0390 0.7308 ± 0.0340 0.3842 ± 0.0178 0.2967 ± 0.0093 0.3445 ± 0.0088 0.3048 ± 0.0080 

ai -0.2946 ± 0.0550 -0.5193 ± 0.0465 -0.4313 ± 0.0432 -0.3921 ± 0.0313 -0.6759 ± 0.0355 -0.6038 ± 0.0339 

as 0.1827 ± 0.0385 -0.0565 ± 0.0374 0.0128 ± 0.0370 0.0129 ± 0.0285 -0.0042 ± 0.0293 -0.0614 ± 0.0294 |-J 

n 31 21 21 30 30 32 

Range 0.8-1.6 0.6-1.8 1.0-3.5 1.3-4.3 1.6-4.8 1.7-5.0 

Median dL (Lq) 0.0180 0.0100 0.0110 0.0089 0.0120 0.0066 



Note. — Tile values , ai, a2, a^, a^, arc tlic coefTieicnts to tlic color-mctallieity-luminosity relations in tire form of Equation^ Tlie bottom tlirec rows 
siiow n, tire number of data points used in the fit, Range, tlie range in color index wlierc eacli relation Irolds true, and Median dL, tlic median absolute deviation 
of the fit. See Figure^for a graphical representation of the solutions, and refer to Section 14.31 for details. 



Table 11. Spectral Type Lookup Table 



Spectral T^pp zt cr zb cr 

Type n (K) (Aq) 



KO 


5 


5347 


± 


20 


0, 


,7956 


± 


0. 


,0076 


Kl 


1 


5147 


± 


14 


0, 


,8051 


± 


0. 


,0035 


K2 


2 


5013 


± 


14 


0, 


.7152 


± 


0, 


,0006 


K3 


2 


4680 


± 


15 


0, 


,7867 


± 


0, 


,0006 


K4 


1 


4507 


± 


58 


0, 


,7390 


± 


0, 


,0190 


K5 


3 


4436 


± 


74 


0, 


,6876 


± 


0, 


,0021 


K7 


3 


3961 


± 


11 


0, 


,6027 


± 


0, 


,0047 


MO 


1 


3907 


± 


35 


0, 


,5773 


± 


0, 


,0131 


MO. 5 


2 


3684 


± 


9 


0, 


,5092 


± 


0, 


,0013 


Ml 


1 


3497 


± 


39 


0, 


,3982 


± 


0, 


,0091 


Ml. 5 


4 


3674 


± 


10 


0, 


,4979 


± 


0, 


,0026 


M2 


1 


3464 


± 


15 


0, 


,3924 


± 


0, 


,0033 


M2.5 


1 


3442 


± 


54 


0, 


,2990 


± 


0, 


,0100 


M3 


3 


3412 


± 


19 


0, 


,4097 


± 


0, 


,0075 


M3.5 


1 


3104 


± 


28 


0, 


,3232 


± 


0, 


,0061 


M4 


1 


3222 


± 


10 


0, 


,1869 


± 


0, 


,0012 


M5.5 


1 


3054 


± 


79 


0, 


,1410 


± 


0, 


,0070 



Note. — 



Sec Section 14. 41 for details. 
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Table 12. Habitable Zone Regions 



Star 


i?i 


nner 




Name 


(AU) 


(AU) 


GJ 15A 


Q 


143 


0.286 


GJ 33 


Q 


454 


0.902 


GJ 53A 


0. 


.583 


1.160 


GJ 75 


0. 


.618 


1.230 


GJ 105 


Q 


471 


0.934 


GJ 144 


Q 


505 


1.003 


GJ 166A 


Q 


559 


1.111 


GJ 205 


Q 


238 


0.475 


GJ 338A 


Q 


252 


0.602 


GJ 338B 


0. 


.243 


0.484 


GJ 380 


0, 


302 


0.601 


GJ 411 


0, 


137 


0.275 


GJ 412A 


0, 


142 


0.284 


GJ 436 





155 


0.311 


GJ 526 


0. 


.183 


0.367 


GJ 551 


0. 


.039 


0.078 


GJ 570A 


0. 


.414 


0.821 


GJ 581 


0. 


103 


0.208 


GJ 631 


0. 


558 


1.111 


GJ 687 


0. 


142 


0.285 


GJ 699 


0. 


057 


0.115 


GJ 702A 


0. 


.623 


1.242 


GJ 702B 


0. 


.359 


0.712 


GJ 725A 


0. 


.121 


0.242 


GJ 725B 


0, 


092 


0.186 


GJ 764 







1.105 


GJ 809 


0. 


215 


0.430 


GJ 820A 


0. 


351 


0.697 


GJ 820B 


0. 


.265 


0.528 


GJ 845 


0. 


.417 


0.828 


GJ 880 


0. 


.218 


0.435 


GJ 887 


0. 


184 


0.368 


GJ 892 


0, 


467 


0.926 


HD 4614 


0. 


907 


1.820 


HD 5015 


l 


507 


3.023 


HD 16895 


1 


194 


2.401 


HD 19373 


1 


207 


2.419 


HD 20630 


0. 


759 


1.519 


HD 22484 




415 


2.838 


HD 30652 


l 


296 


2.620 


HD 34411 


l 


[)92 


2.185 


HD 39587 


0. 


846 


1.697 


HD 48737 


2, 


633 


5.323 


HD 56537 




532 


7.340 


HD 58946 


1. 


.748 


3.568 


HD 81937 


2, 


.943 


5.971 


HD 82328 


2, 


210 


4.454 


HD 82885 


0. 


.756 


1.507 


HD 86728 


0. 


.987 


1.970 


HD 90839 


1. 


.005 


2.023 


HD 95418 


4. 


603 


9.811 


HD 97603 


3. 


291 


6.861 


HD 101501 


0, 


676 


1.345 


HD 102870 


1. 


513 


3.042 


HD 103095 


0. 


.416 


0.825 


HD 109358 


0. 


.898 


1.796 


HD 114710 


0. 


.950 


1.906 


HD 118098 


2. 


.740 


6.722 


HD 126660 


1. 


.607 


3.236 


HD 128167 


1, 


.424 


2.883 


HD 131156 


0. 


677 


1.351 


HD 141795 


2. 


294 


4.776 


HD 142860 


1. 


.374 


2.769 


HD 146233 


0. 


.879 


1.761 


HD 162003 


2. 


039 


4.093 


HD 164259 


1. 


.926 


3.896 



- 54 - 



Table 12 — Continued 



star 


dinner 


Houtor 


Name 


(AU) 


(AU) 


HD 173667 


1.962 


3.959 


HD 177724 


3.658 


7.772 


HD 182572 


1.141 


2.284 


HD 185395 


1.614 


3.257 


HD 210418 


3.307 


6.867 


HD 213558 


3.173 


6.733 


HD 215648 


1.734 


3.488 


HD 222368 


1.487 


2.997 



Note. — Habitable zone bound- 
aries for the K- and M-stars in this 
work and A- , F- , and G-st ars from 
iBoyajian et all i2012|) (top and bot- 
tom portion of the table, respec- 
tively). See Section[7]for details. 
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Fig. 1. — Calibrated observations plotted with the limb-darkened angular diameter fit for 
each star. 
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Fig. 2. — Calibrated observations plotted with the limb-darkened angular diameter fit for 
each star. 
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Fig. 3. — Calibrated observations plotted with the limb-darkened angular diameter fit for 
each star. 



-58- 



1.0 

0.8 
I 0.6 
> 0.4 

0.2 

0.0 

M 0.2 

§ 0.1 

^ 0.0 

S -0.1 

DC -0.2 







r GJ809 












f i 






140 


160 180 
Baseline/Wavelength X 1 0"^ 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 


200 



^ -0.1 
-0.2 




100 120 140 180 180 200 
Baseline/Wavelength X 1 0"^ 




110 120 130 140 150 
Baseline/Wavelength X 1 0^^ 



160 



Fig. 4. — Calibrated observations plotted with the limb-darkened angular diameter fit for 
each star. 
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9 Berger et al. 2006 

van Belle & von Broun 2009 

Lone et al. 2001 

Demory et al. 2009 



1.0 1.2 

(this work) 



1.4 



1.6 



Limb- darkened diameters of the stars observed in common with the works 



of iLane et ahl (120011): ISegransan et al.l ( 120031 ): iBerger et ahl (120061 ): Ivan Belle fc von Braun 
(120091 ): iDemory et al.l ( l2009l ) compared to the values measured in this work (left to right: 
GJ 699, GJ 15A, GJ 380, GJ 411, GJ 166A). Results of this consistency check leads to an 
agreement of 6'LD,this work/6'LD, reference = 1-008. See Section [2] for details. 
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Fig. 6. — Empirical color- metallicity-temperature relations presented in Section 14.11 The 
color of the data point reflects the metallicity of the star as depicted in the legend. The col- 
ored lines are solutions to the metallicity dependent fits, where the line color (orange, green, 
teal) represents our solution for an iso- metallicity liii e to [F e/H] = 0.0, —0.25, —0.5. The 
dash-dotted lines are the solutions from lLejeune et al. (|l998h . Solut ions from the empirical 
relation established for dwarfs via interferometric measurements in van Belle fc von Braun 



(120091 ) is shown as a dashed-triple-dotted line. Color-temperature curves from the BT-Settl 
PHOENIX model atmospheres at a [Fe/H] = are illustrated as a dotted line. See Equa- 
tion [3] for the form of the equation and Table M for the coefficients and statistical overview 
for each color-metallicity-temperature solution in this work. 
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Fig. 7. — Temperatures determined in this work compared to those found in the hterature 
for each star. The position in the x-direction is our measured temperature and associated 
1-cr error on temperature. The extent of the y-error bars indicate the range in temperature 
that can be found in the hterature for each star. See Section 14.11 for details. 
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Fig. 8. — Empirical color-metallicity-radius relations presented in Section 14.21 The color 
of the data point reflects the metallicity of the star. The colored lines are solutions to 
the metallicity dependent fits, where the line color (red, orange, green, teal) represents our 
solution for an iso-metallicity line to [Fe/H] = +0.25, 0.0, —0.25, —0.5. The panels displaying 
the (B — V), (V — J), (V — H), and (V — K) relations also include Dartmouth 5 Gyr isochrones 
as dash-dotted ([Fe/H] = 0.0) and dotted lines ([Fe/H] = -0.5). Within the {B - V) panel, 
we include solutions with the same Dartmouth model specifications, but using the semi- 
empirical color transformation (blue lines) to derive the [B — V) colors. See Equation H] for 
the form of the equation and Table M for the coefficients and statistical overview for each 
color-metallicity-radius solution. 
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Fig. 9. — Empirical color- metallicity-luminosity relations presented in Section I^l3l The color 
of the data point reflects the metallicity of the star. The colored lines are the solutions for 
the metallicity dependent fits, where the line color (red, orange, green, teal) represents our 
solution for an iso-metallicity line to [Fe/H] = +0.25, 0.0, —0.25, —0.5. The panels displaying 
the (B — V), (V — J), (V — H), and (V — K) relations also include Dartmouth 5 Gyr isochrones 
as dash-dotted ([Fe/H] = 0.0) and dotted lines ([Fe/H] = -0.5). Within the {B - V) panel, 
we include solutions with the same Dartmouth model specifications, but using the semi- 
empirical color transformation (blue lines) to derive the [B — V) colors. See Equation [5] for 
the form of the equation and Table [10] for the coefficients and statistical overview for each 
color-metallicity-luminosity solution. 
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Fig. 10. — The 5 Gyr models isochrones of Padova, Dartmouth, BCAH98, and Y onsei- Yale 
Jcirardi et al.lboooi IPotter et allboosi : ISaraffe et aPllQQsl : IPemarque et al.ll2004l ) displayed 
along with the single star measurements presented in Table [6] (filled circles). In the top-right, 
top- left, and bottom- left plots, the dashed line corresponds to the solar metallicity model 
(see color legend within plots), and the hashed region spans a metallicity of to —0.5 dex. 
In the bottom-right plot, we show observations of single stars (filled circles; Table E]) and 
binary stars (open squares, see Section l57^ . The solid and dashed lines mark a change in 
metallicity from [Fe/H] = to —0.5. Characteristic one-sigma mass errors for single stars 
[aMss) are shown close to the bottom axis. See Section |5] for details. 
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Fig. 11. — Top: Temperature-luminosity relation derived here is shown as a solid black line 
(Equation |6]) for the approximate range the relation holds true. Dartmouth model 5 Gyr 
isochrones are also shown as dash-dotted ([Fe/H] = 0) and dotted lines ([Fe/H] = —0.5). 
Middle: The fractional difference in luminosity (observed versus our fit to the data) is shown 
versus temperature. The dotted line represents a zero deviation from the fit. In both the top 
and middle plots, errors are not plotted, but typically smaller than the size of the data points. 
Bottom: The observed temperature versus the temperature predicted by interpolating the 
observed luminosity in the custom-tailored Dartmouth model for each star (Tmoci)- The 
dotted line shows a 1:1 relation, and the y-errors shown are the scaled observational errors 
in temperature. In all plots, the color of the data point reflects the metallicity of the star as 
depicted in the legend. See Section 15.11 for details. 
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Fig. 12. — Top: Radius-luminosity relation derived here is shown as a solid black line 
(Equation [7]) for the approximate range the relation holds true. Dartmouth model 5 Gyr 
isochrones are also shown as dash-dotted ([Fe/H] = 0) and dotted lines ([Fe/H] = —0.5). 
Middle: The fractional difference in luminosity (observed versus fit to the data) is shown 
versus radius. The dotted line represents a zero deviation from the fit. In both the top and 
middle plots, errors are not plotted, but typically smaller than the size of the data points. 
Bottom: The observed radius versus the radius predicted by interpolating the observed 
luminosity in the custom-tailored Dartmouth model for each star (i?Mod)- The dotted line 
shows a 1:1 relation, and the y-errors shown are the scaled observational errors in radius. 
In all plots, the color of the data point reflects the metallicity of the star as depicted in the 
legend. See Section 15.21 for details. 
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Fig. 13. — Temperature-radius relation is shown as a solid black line (Equation |H]) for the 
approximate range that the relation holds true. The relation including the Sun (Equation [H]) 
is shown as a blue line, where the dotted line is an extrapolation of the curve to higher 
temperatures with the Sun as a reference point. The color of the data point reflects the 
metallicity of the star. Errors are not plotted, but are typically of the sizes of the data 
points. The left-most orange dot is the Sun. Solar metallicity Dartmouth models are used 
to plot the positions at 1 Gyr for stars with masses of 1.0 M©, O.QM©, and O.SM© (black x, 
red +, and blue * respectively). Arrows point to their positions at an age of 4.5 Gyr. See 
Section [5l3] for details. 
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Fig. 14. — Top: The observed temperature versus the temperature predicted by interpolating 
the observed radius in the custom-tailored Dartmouth model for each star (TMod)- Bottom: 
The observed radius versus the radius predicted by interpolating the observed temperature 
in the custom-tailored Dartmouth model for each star (-Rmoci)- The dotted line shows zero 
deviation from the model, and the y-errors shown are the scaled observational errors in 
temperature (top plot) and radius (bottom plot). In all plots, the color of the data point 
reflects the metallicity of the star as depicted in the legend. See Section 15.31 for details. 
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Fig. 15. — Mass-radius relations for single and binary stars as expressed in Equations [10] 
and [TTJ The filled circles and solid line are the data and solution for single stars. The open 
squares and dotted line are for the EB stars. The measured 1-sigma errors are shown for 
radii, but are typically smaller than the data point. Single star mass errors are not shown 
for clarity. We show a typical single star mass error bar for a given mass at the bottom of 
the plot window indicating a value of aM^s ~ 10 %. Although the mass errors for single 
stars are large, we do not detect any metallicity dependence on the mass-radius relation. See 
Section 15.41 for details. 
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Fig. 16. — Fractional offset in radii versus mass {top), metallicity [Fe/H] [middle) and nor- 
malized X-ray to bolometric luminosity Lx/Lbol {bottom) for single stars compared to 5 Gyr 
Dartmouth model isochrones custom-tailored to each star. Characteristic one-sigma mass 
errors for single stars {aMss) are shown close to the bottom axis in the top plot. The extent 
of the y-errors in all plots show a fractional deviation in radius of 10% (due to errors in mass 
estimates). The left middle and bottom panels show the full sample of K- and M-dwarfs 
(M < 0.9 Mq). The right middle and bottom panels show only the M-dwarfs (M < 0.6 Mq), 
where the blue points indicate stars with masses < 0.42 Mq. The dotted line indicates a 
zero deviation. See Section [5751 for details. 
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Fig. 17. — Stellar radius versus temperature for single stars (red circles) and EBs (blue 
squares). The solid line is the polynomial fit to the empirically determined single star data 
presented in Equation [91 The solar metallicity solution from the Dartmouth models is shown 
as the dash-dot line. The tabulated values for temperatures and radii for main sequence stars 
in Allen's AQ (jCoxll2000l ) is plotted as a dashed line, predicting a decreased temperature of 
~200 to 300 K for a star of given radius. Note the correlati ons to binar y star temperatures 
and radii show to be much more suited for the Allen's AQ ICoxl (120001 ) temperature scale. 
See Section [6] for details. 
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Fig. 18. — An empirically determined H-R diag ram for the collection of stars presented in 
this paper and in Paper I (IBoyajian et al.ll2012l ). The yellow five-pointed star is the Sun. 



